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Abstract 



o 

' We survey some topics involving the Whitham equations, concentrating on the role of iptp* 

(or square eigenfunctions) in averaging. Some connections to symplectic geometry and Seiberg- 
Witten theory are indicated. 

(N ; 

1 INTRODUCTION 

. A preliminary sketch of some of this material was given in ]15| and in view of several recent devel- 

opments (cf. |, [fj, m H, H H> HI. H H E§ H) we decided t0 expand and rewrite [[5) in 



two parts, of which this is the first. The second part will deal with theory of kernels related to ■0V'* 

There are many physical situations where averaging procedures of various types are employed 
and we refer to @ for a general background (see || fH |§ H 11 M H II 11 H H' 11 11 
for weakly deformed soliton lattices, multiphase averaging, etc.). The Whitham equations arise in 
this context and we will look at some of this here; we attempt to organize a subset of information 
in a coherent and more or less rigorous manner, while clarifying various formulas and assertions in 
the literature. Connections to dispersionless theory are indicated and material has been selected 
which will be useful in further studies of topological field theory (TFT), Landau-Ginzburg (LG) 
theory, Seiberg-Witten (SW) theory, etc. The exposition in this first paper is mainly expository 
but we indicate some new connections and results (cf. also |L6|, [l?], [18], |2(J). The principal theme 
involves averaging of square eigenfunctions for KdV or averaging of quantities related to ipip* for 
KP, where ip is the Baker- Akhiezer (BA) function; this point of view (along with general use of the 
BA function) is emphasized for derivation of Whitham equations, construction of differentials, and 
moduli considerations. Many basic formulas involving theta functions, differentials, etc. are also 
exhibited in order to make the text useful as a launching pad for further investigation. 

2 CURVES AND MODULI 
2.1 Riemann surfaces and BA functions 

Following 0, H |[ [H H 11, i§[ |]| |5£| Q we take an arbitrary Riemann surface S of genus g, 
pick a point Q and a local variable 1/k near Q such that k(Q) = 00, and, for illustration, take 
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q(k) = kx + k 2 y + k 3 t. Let D = Pi + ■ ■ ■ + P g be a non-special divisor of degree g and write ip for 
the (unique up to a constant multiplier by virtue of the Riemann-Roch theorem) Baker-Akhiezer 
(BA) function characterized by the properties (A) ip is meromorphic on £ except for Q where 
ip{P)exp{-q{k)) is analytic and (*) ip ~ exp(q(k))[l + J2T(^/ kj )} near Q- ( B ) 0n S /<3> ^ has 
only a finite number of poles (at the Pi). In fact ip can be taken in the form (P € S, Pq ^ Q) 



„, i, P) = e,p[ f (^ + ^ + *«»)] • 9( " 4(P) o^P^ + W + Z ° ) ^ 

Jpo e(^(P) + z ) 

where dfi 1 = dfc + • • • , df7 2 = d(fc 2 ) + • • • , dft 3 = d(k 3 ) + ■ ■■ , Uj = J B , dft\ Vj = J B .dfl 2 , Wj = 
j B dVl 3 (j = l,---,g), zq = —A(D) — K, and O is the Riemann theta function. The symbol 
~ will be used generally to mean "corresponds to" or "is associated with"; occasionally it also 
denotes asymptotic behavior; this should be clear from the context. Here the dVtj are meromorphic 
differentials of second kind normalized via J A dflj — (A,, Bj are canonical homology cycles) and 

we note that xdVL 1 +ydQ 2 +tdQ 3 ~ dq{k) normalized. A ~ Abel-Jacobi map A{P) — (Jp o du>f.) where 
the dcuk are normalized holomorphic differentials, k = 1,- • ■ ,g, J A duik — djk, and K = (Kj) ~ 
Riemann constants (2K = —A(Ks) where Ks is the canonical class of E ~ equivalence class of 
meromorphic differentials). Thus Q(A(P) + zq) has exactly g zeros (or vanishes identically. The 
paths of integration are to be the same in computing f p dVf 1 or A(P) and it is shown in ||, ||, |l5|, |2^| 
that ip is well defined (i.e. path independent). Then the £j in (*) can be computed formally 
and one determines Lax operators L and A such that d v ip = Lip with dtip = Aip. Indeed, given 
the £j write u = —2d x £_i with w — 3£,id x !;i — 39 2 ^i — 3d x £_2- Then formally, near Q, one has 
(-d y +dl + u)ip = 0(l/k)exp(q) and {-d t +d 3 + (3/2)ud x + w)ip = 0(l/k)exp(q) (i.e. this choice 
of u, w makes the coefficients of k n exp(q) vanish for n = 0,1,2,3). Now define L = d 2 + u and 
A = d 3 + (3/2)ud x + w so d y ip — Lip and dtip = Aip. This follows from the uniqueness of BA 
functions with the same essential singularity and pole divisors (Riemann-Roch). Then we have, via 
compatibility L t — A y — [A, L], a KP equation (3/4)u yy = d x [u t — {1/ A){&uu x + u xxx )\ and therefore 
such KP equations are parametrized by nonspecial divisors or equivalently by points in general 
position on the Jacobian variety </(£)■ The flow variables x, y, t are put in by hand in (A) via q(k) 
and then miraculously reappear in the theta function via xll + yV + tW; thus the Riemann surface 
itself contributes to establish these as linear flow variables on the Jacobian. The pole positions Pi 
do not vary with x, y, t and (f ) u — 2d 2 logQ{xU + yV + tW + zo) + c exhibits 8 as a tau function. 

We recall also that a divisor D* of degree g is dual to D (relative to Q) if D+D* is the null divisor 
of a meromorphic differential dCl = dk + (f3/k 2 )dk + ■ ■ ■ with a double pole at Q (look at £ = 1/k to 
recognize the double pole). Thus D + D* - 2Q - so A{D*)-A{Q)+K = ~[A{D)- A(Q)+K}. 
One can define then a function ip* (x,y,t, P) = exp(—kx — k 2 y — k 3 t)[l + £*/fc) + • • •] based on 
D* (dual BA function) and a differential dCl with zero divisor D + D* , such that (p = ipip*dti is 
meromorphic, having for poles only a double pole at Q (the zeros of dQ cancel the poles of ipip*). 
Thus ipip*dCl ~ ipip*(l + (P/k 2 + ■ ■ -)dk is meromorphic with a second order pole at oo, and no other 
poles. For L* = L and A* = -A + 2w- (3/2)u x one has then (d y + L*)ip* = and (d t + A*)ip* = 0. 
Note that the prescription above seems to specify for ip* (U = xU + yV + tW, Zq = — A(D*) — K) 

- (x<m 1 +ydQ 2 +tdn 3 ) ®(A(P) -U + z*) 

In any event the message here is that for any Riemann surface E one can produce a BA function 
ip with assigned flow variables x, y, t, • • • and this ip gives rise to a (nonlinear) KP equation with 
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solution u linearized on the Jacobian J(S). For averaging with KP (cf. f«], we can use 
formulas (cf. ([H]) and (gg)) 



tp = e px+Ev+nt ■ cj)(Ux + Vy + Wt, P) (2.3) 

^* = e - P x-Ey-nt . ^^_ Ux ~Vy- Wt, P) (2.4) 

to isolate the quantities of interest in averaging (here p — p(P), E = E(P), ft = £l(P), etc.) We 
think here of a general Riemann surface S ff with holomorphic differentials diok and quasi-momenta 
and quasi-energies of the form dp = dO, 1 , dE = dQ 2 , dfl — dfl 3 , ■ ■ ■ (p — f p ^ dfl 1 etc.) where the 
dQi = dVlj = d(X' + 0(X^ 1 )) are meromorphic differentials of the second kind. Following one 
could normalize now via 3? J A d£l k = dQ = 0. Then write e.g. Uk = (l/2ni) § A ^ dp and 

Uk+g — — (l/2m) § Bk dp (k — 1. • • • ,g) with similar stipulations for Vk ~ § dfl 2 , Wk ~ § dfl 3 , etc. 
This leads to real 2g period vectors and evidently one could also normalize via § A dVl k — or 
^ § A dVl k =^s§ B dfl k = (further we set Bjk = § Bk duij). 



2.2 General remarks on averaging 



Averaging can be rather mysterious at first due to some hasty treatments and bad choices of notation 
plus many inherent difficulties. Some of the clearest exposition seems to be in |37], [||, ^36], |7g| 



whereas the more extensive developments in e.g. 
until one realizes what is going on. 
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| ||^, Q may be confusing at first, 



One choice of background situation involves an examination of dispersionless limits with no 
concern for periodicity or Riemann surfaces (cf. here j^, O, [l3|, [l4|, |4lJ [I?], f74|). Thus given e.g. a 



KdV equation u t + 6uu x 



= set ex = X and et = T leading to u|, + 6u £ it ( 



- + e 2 u 



xxx 



0. The 



study of u e — > u where u satisfies the Euler equation ut + 6uux = is a very delicate matter, of great 
interest in applications and in PDE but we do not discuss this here (cf. [pOl |7(| ) . On the other hand 
the purely algebraic passage of the background mathematics of KdV (involving Lax operators, the 
KdV hierarchy, tau functions, vertex operators, etc.) to the corresponding background mathematics 
of the dispersionless theory, is relatively easy and will be indicated below. Moreover it is of great 
importance in an entirely different direction, namely in the study of topological field theory, strings, 
and 2-D gravity (see e.g. §, |, H |J, ^ || lH HI HI ll HI ) ■ We will insert such 
material later as appropriate. 

Let us follow |33| (cf. also (5^]) in order to have a suitably complicated example leading directly 
to matters of interest here, so as background we consider the KP equation 



-a 2 u y 



d x [u t 
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(2.5) 



(a 2 = 1 - KP-2; a 2 = -1 - KP-1). Here crL - d 2 - u with ad y i/j = (d 2 x - u)tp and A ~ 
d 3 — (3/2)ud x + w, which is slightly different than before (to connect notations let u — > — u and insert 
a). Then one has the compatibility equations [d y — L, dt — A] = as before and it should be noted 
that for any function g(t), t ~ x,y,t,ti,- ■ •, operators L — gLg~ 1 +d y g g^ 1 and A — gAg~ 1 +dtg g 
could be used, corresponding to a new wave function tp = gip. The notation of |53j also involves 



differentials d&h 



dk l (l + 0(1/ k)] and J A dtlj = with U k = (1/2™) J Bk dflj 
in (^) with potentials u given via (f). We recall that Riemann theta functions are more precisely 



dVt 1 where d£li 

Tk _ 



(so U k ~ Uk, U$ ~ Vk, and U£ ~ Wk, up to factors of 2iri). There will be solutions ip as before 
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written as Q(z\B) where B ~ a certain matrix (cf. below), and z ~ (zi, • • • , z g ) ~ xU + yV + tW + zq 
for example. One knows via || |30| that u of the form (f) u — 2d^logQ(z\B) + c is a solution of 
KP if and only if the matrix B defining © is the B-period matrix of a Riemann surface determined 
via f B duij (with z as defined). Now consider the spectral theory of 2-D periodic operators (**) 
(crdy — d 2 + u(x, y))ip = where u(x, y) — u[x + a±, y) — u[x, y + 02). Bloch solutions are defined 
via 

ip(x + ai,y, wi,w 2 ) = Wiip(x, y, W!,w 2 ); (2.6) 

ip(x, y + aa,wi, w 2 ) = w 2 ip(x, y, w u w 2 ) 

and we assume ^(0, 0, W\, W 2 ) = 1. The pairs Q = (w\,w 2 ) for which there exists such solutions is 
called the Floquet set T and the multivalued functions p(Q) and E(Q) such that w\ = exp(ipai) 
with w 2 = exp(iEa 2 ) are called quasi-momentum and quasi-energy respectively. By a gauge trans- 
formation ip — > exp(h(y))ip, with d y h(y) periodic one obtains solutions of (**) with a new potential 
u = u — adyh so we can assume J Q 1 u(x, y)dx = 0. 

For Mq = ad y — d 2 (with u = 0) the Floquet set is parametrized by k G C such that 
Wi = exp{ika\) and w® — exp(—k 2 a 2 /a) and the Bloch solutions are ^{x, y, k) — exp(ikx — k 2 y/a). 
The adjoint Bloch solutions are ip + (x, y, k) — exp(—ikx + k 2 y/a) satisfying (ad v + d 2 )ip + = 0. The 
image of the map k — > (u^, w%) £ C 2 is the Floquet set for Mo corresponding to the Riemann surface 
with intersections corresponding to pairs k 7^ k' such that w®(k) = w^(k'), i = 1,2. This means 
k — k' = (2irN/ai) and k 2 ~(k') 2 — (2iriaM / a 2 ) where N, M are integers, creating "resonant" points 
k = fcjv,M = i^N/ai)— (iaMa±/Na 2 ), N 0, k' = k-N.-M- Then for kg ^ fcjv,M and u sufficiently 
small one can construct a formal Bloch solution of (**) in the form of a convergent perturbation 
series (for any a). Thus outside of some neighborhoods of the resonant points one can obtain a 
Bloch solution tp(x,y,ko) which is analytic in fco, but the extension of tp to a resonant domain can 
be tricky. For 3?ct = (as in KP-1) the resonant points are dense on the real axis whereas for 5Rcr 7^ 
(as in KP-2) there are only a finite number of resonant points in any finite domain of C. In the latter 
case one can glue handles between domains around resonant points and create a Riemann surface 
r of Bloch solutions ip(x,y,Q), Q £ T. Moreover, if the potential u can be analytically extended 
into a domain \Qx\ < t\, < t 2 , then the technique can also be adapted even when u is not 
small. Such Bloch solutions, normalized by -0(0, 0,Q) = 1, are meromorphic on T and in the case 
of a finite number of handles a one point compactification of T is obtained so that ip is in fact the 



BA function for T. Generally speaking finite zone situations as in (2.1) with potentials given by 
Riemann theta functions are quasi-periodic in nature. To single out conditions for periodicity one 
asks for meromorphic differentials dp and dE on Y having their only singularities at Q ~ point at 
00 of the form dp — dk(l + 0(k~ 2 )) and dE = ia~ 1 dk 2 (l + 0(k~ 3 )), normalized so that all periods 
are real, and satisfying, for any cycle C on T, § c dp = (2imc / ai) with § c dE = {2-Kmc / a 2 ) where 
n Ci m c are integers. Then the corresponding potentials u(x,y) will have periods ai and a 2 in x 
and y, with multipliers wi(P) = exp(ia \ f dp) and w 2 (P) = exp(ia 2 J dE). We go now to finite 
zone (or quasiperiodic) situations as in ( |2.l| ) with potentials as in (f), where z ~ xlJ + yV + tW + z 
can be written as (zk) — [x J b ft 1 + y J b Q 2 + t J, Q 3 + z%\ — (Cfc + ^ ~ ^-1'" ->9- Since 
Q(z + 2niN) — O(z) we could set ( = id so that as a function of 9 = (9k), u is periodic of period 
2ir in each variable 0k ■ Now one wants to consider modulated finite zone situations where solutions 
are of the form u = u (xU + yV + tW\I) — Uo(#i, ■ ■ • , 9 g \I\, ■ ■ ■ , I n ) where uq is periodic in the 9j 
with U, V, W = U, V, W(I). One assumes there will be slow variables X = ex, Y = ey, and T = et 
with fast variables x,y,t so that "asymptotic" solutions 

u = u (^S(X, Y, T)\I(X, F, T)) + eu x {x, y, t) + e 2 u 2 (x, y,t) + --- (2.7) 
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can be envisioned. In practice the parameters Ik depend on the moduli of our Ricmann surface, which 
are then allowed to change smoothly with the slow variables X, Y, T (so U, V, W also can change), 
and one looks for solutions ( |2.7| ) with uniformly bounded u\. We do not look for convergence of the 
series in (2.7) nor check any other features of "asymptotic solution" or "asymptotic series". Such 
procedures are standard in the study of what are called weakly deformed soliton lattices and we 
emphasize that the introduction of slow variables here is an assumption based on modulated wave 
trains. 



2.3 Hyperelliptic curves and KdV averaging 

One knows that hyperelliptic curves play a special role in the theory of algebraic curves and Riemann 
surfaces (see e.g. [§§|[|§^],|§|§||,|§,||,H,|9|,^|). For hyperelliptic Riemann surfaces 
one can pick any 2g + 2 points Xj G P 1 and there will be a unique hyperelliptic curve E g with a 
2-fold map / : E g — > P 1 having branch locus B = {Xj}. Since any 3 points Ai, Xj, Xk can be sent to 
0, 1, oo by an automorphism of P 1 the general hyperelliptic surface of genus g can be described by 
(2<7 + 2) — 3 = 2g — 1 points on P 1 . Since / is unique up to an automorphism of P 1 any hyperelliptic 
S ff corresponds to only finitely many such collections of 2g — 1 points so locally there are 2g — 1 
(moduli) parameters. Since the moduli space of algebraic curves has dimension 3g — 3 one sees that 
for g > 3 the generic Riemann surface is nonhyperelliptic whereas for g = 2 all Riemann surfaces 
are hyperelliptic (with 3 moduli). For g = 1 we have tori or elliptic curves with one modulus r and 
g = corresponds to P 1 . In many papers on soliton mathematics and integrable systems one takes 
real distinct branch points Xj, 1 < j < 2g + 1, and oo, with Ai < A2 < • ■ ■ < A2 9 +i < 00 and 

23+1 

M 2 = ]J(X-X j ) = P 2g+1 (X,X j ) (2.8) 
1 

as the defining equation for S g . Evidently one could choose Ai = 0, A2 = 1 in addition so for 3 = 1 
we could use < 1 < u < 00 for a familiar parametrization with elliptic integrals, etc. One can take 
dX/ ix, XdX/fi, • • • , A 9_1 rfA//i as a basis of holomorphic differentials on T, g but usually one takes linear 
combinations of these denoted by du)j, 1 < j < g, normalized via § A dujj = Sij, with period matrix 
defined via § B du>j = Ely. The matrix II = (II jj) is symmetric with 311 > and it determines 
the curve. Frequently in situations arising from KdV (Korteweg-deVries equation) for example one 
regards the intervals [Ai, A2], • • • , [A2g+i, 00) as spectral bands and intervals (A2, A3), • • • , (X-2 g , A2 9 +i) 
as gaps with the ai cycles around (A2;, A2i+i) (i — 1, • • • ,g). One will also want to consider another 
representation of hyperelliptic curves of genus g via 

2g+l 

M 2 = H(X-X J ) = P 2g+2 (X,X J ) (2.9) 


where 00 is now not a branch point and in fact there are two points \x± corresponding to A = 00. 

We recall next some of the results and techniques of |37], |66t w here one can see explicitly the 
nature of things. The presentation here follows First from |66|, in a slightly different notation, 
write q t = 6qq x — <lxxx with Lax pair L — —d^ + q, B = -Adl + 3{qd x + d x q), L t = [B,L], LtJj — Xtp, 
and ipt = Bip. Let tj) and <fi be two solutions of the Lax pair equations and set \& = ip4>] these are 
the very important "square eigenfunctions" which arise in many ways with interesting and varied 
meanings (cf. 0, |D], |TJ, Q). Evidently * satisfies 

[-flg + 2{qd x + d a q)]V = 4Xd x V; ft* = -2q x ^> + 2(q + 2A)9 X * (2.10) 
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From ( |2 . 1 OP one finds imm ediat ely the conservation law (C): d t [V] + d x [6(q - 2A)* - 2d^l] = 0. If 
one looks for solutions of ( [2.1C| ) of the form ^>(x,t, A) = 1 + Y^Ti^j( x ^)]^~^ as A ^ oo then one 
obtains a recursion relation for polynomial densities 

d^i+i = \-\dl + (qd x + d x q)}*j (j = 1, 2, • • •); *o = 1 (2.11) 

Now consider the operator L = —d% + q in L 2 (—oo, oo) with spectrum consisting of closed intervals 
separated by exactly N gaps in the spectrum. The 2N + 1 endpoints A^ of these spectral bands 
are denoted by — oo < Ao < Ai < • • • < X 2 n < oo and called the simple spectrum of L. They 
can be viewed as constants of motion for KdV when L has this form. We are dealing here with 
the hyperelliptic Riemann surface determined via -R 2 (A) = Hq^A — ^fc) (°° 1S a branch point) and 
one can think of a manifold Ai of TV-phase waves with fixed simple spectrum as an TV-torus based 
on 9j £ [0, 2tt). Hamiltonians in the KdV hierarchy generate flows on this torus and one writes 
q = <7at(#i, • • • , 9n) (the 9 variables originate as in our previous discussion if we use theta functions 
for the integration - cf. also below). Now there is no y variable so let us write 9j = xKj +twj (we will 
continue to use du>j for normalized holomorphic differentials) For details concerning the Riemann 
surface we refer to § |, ^, || and will summarize here as follows. For any q^ as indicated 
one can find functions /ij(x,t) via fy(x,t, A) = Yii ~ f l j( x ^)) where /ij(x,t) £ [X^j-i, X^j] and 
satisfies 

d xf Mj = -2i(R(fMj)/ - Mi)); (2.12) 

2N 

dtfij = -2i[2(J2 Xk-2^2 fa)] ■ (R(fXj)/ JJ( M j - fa) 

ijij i=£j 

In fact the fij live on the Riemann surface of R(X) in the spectral gaps and as x increases i^j travels 
from A2J-1 to X%j on one sheet and then returns to X 2 j-i on the other sheet; this path will be called 
the j th /i-cycle (~ Aj). In the present context we will write the theta function used for integration 
purposes as 0(z, r) = J2mez N exp[ni(2(m, z) + (m, nil))] where z £ Z N and r denotes the N x N 
period matrix (r is symmetric with 3r > 0). We take canonical cuts Ai, Bi (i = 1, • • • ,N) and 
let duij be holomorphic diffentials normalized via J a dio^ = Sjk (the cycle Aj corresponds to a loop 
around the cut Aj). Then q^ can be represented in the form 



2N 

q N {x,t) = A + T -2dlloge{z(x,t);T); A = ^Aj-; (2.13) 



( T ij) = (f du 3 ); r*j = -nj] r = -2V* Xdujj 

JB, i J Aa 



andz(x,t) = — 2i[c N (x— xo)+2(Ac Ar +2c Ar ^tj+d where (c N )i = CiM arises from the representation 
duji = Q2i <kj)j~ l )[dX/ R{X)] (d is a constant whose value is not important here). Then the wave 
number and frequency vectors can be defined via k — —Uttt^c 1 ^ and w — — 8iirr~ 1 [Ac N + 2c N ~ 1 ] 
with 9j(x, t) — KjX + Wjt + 9j (where the 9® represent initial phases). 

To model the modulated wave now one writes now q = qN(9i, • ■ ■ , 9n; A) where Xj ~ Xj {X, T) and 
A ~ {Xj). Then consider the first 2iV+l polynomial conservation laws arising from ( 2.1 0| ) - ( 2.11 ) and 
C for example (cf. below for KP) and write these as d t Tj(q)+d x Xj(q) = (explicit formulas are given 



G 



roughly as follows). We note that the adjoint linear KdV equation (governing the evolution 



of conserved densities) is dtjj + d x ^fj — 6qd x ^j = (7^ ~ Vifj) and (2.11) has the form djj+i = 
(— (l/2)<9 3 + qd + dq)jj. One then rewrites this to show that 6qd x jj — d x [6jj + i — €>q^j + 3d 2 jj] so 
that the adjoint equation becomes 

d tlj + + 6 qij - 6 7i+1 ] = (2.14) 



which leads to ( ^.15| ) and (2.19) below (after simplification of (2.14)). Then for the averaging step 



write d t = edx, etc., and average over the fast variable x to obtain 

d T < T 3 {q N ) > +d x < Xj{q N ) >= (2.15) 
(( p. 15 ) makes the first order term in e vanish). The procedure involves averages 



1 " L 



< Tj(q N ) >= lirriL-,00— / Tj(q N )dx (2-16) 

<J — L 

for example (with a similar expression for < Xj(q^) >) and an argument based on ergodicity is 
used. Thus if the wave numbers Kj are incommensurate the trajectory {q^(x,t); x G (— 00,00)} 
will densely cover the torus M.. Hence we can replace x averages with 

<T M>= J^jNj ■■■J r ^{6))J{de j (2.17) 

For computational purposes one can change the 9 integrals to fj, integrals and obtain simpler calcu- 
lations. By this procedure one obtains a system of 2N + 1 first order partial differential equations 
for the 2N + 1 points Xj(X,T), or equivalently for the physical characteristics (k(X,T),w(X,T)) 
(plus < qjy >). One can think of freezing the slow variables in the averaging and it is assumed that 
( [2.15 ) is the correct first order description of the modulated wave. 



The above argument may or may not have sounded convincing but it was in any case rather 
loose. Let us be more precise following J37j. One looks at the KdV Hamiltonians beginning with 
H = H(q) = iimi_ >00 (l/2£) J^ L (q 2 + (l/2)q 2 )dx (this form is appropriate for quasi-periodic situa- 
tions). Then q t = {q,H} where {/,<?} = Hmi_ KX3 (l/2i) j\(5f /5q)d x (8g/6q)dx (averaged Gardner 
bracket). The other Hamiltonians are found via 

d S J^±l = ( qd + dq - >)^i (m > 0); ^ = 1 (2.18) 
oq 2 oq oq 

where 7j ~ VHj ~ (SHj/Sq) (cf. here |], [3?]). It is a general situation in the study of symmetries 
and conserved gradients (cf. ||10[ ) that symmetries will satisfy the linearized KdV equation (dt — 
6d x q + d^)Q — and conserved gradients will satisfy the adjoint linearized KdV equation (dt — 
6qd x + d x )Q^ = 0; the important thing to notice here is that one is linearizing about a solution q of 
KdV. Thus in our averaging processes the function q, presumed known, is inserted in the integrals. 
This leads then to 



with ( [2.15 ) holding, where < d 2 4> >= implies 



i A / / a / / 

<X j >=lim L ^ 00 — (-6-J±± + Qq N -J-)dx (2.20) 
2L J_ L dq N dq N 
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In p7[ | the integrals are then simplified in terms of \i integrals and expressed in terms of abelian 
differentials. This is a beautiful and important procedure linking the averaging process to the 
Riemann surface and is summarized in Ipral as follows. One defines differentials 



= --A 



'R(X) 

T 



2 J; V 

where the Cj, Ej are determined via $ b . &i = = § h ,&2 (* = 1, 2, • ■ • , N). Then it can be shown 
that 

< * >~< r >~ V < A" >~ V (2.22) 

with f2j ~< T > (d£/£ 2 ) and < A" > (cZ^/C 2 ) ~ 12[(d£/£ 4 ) - 2 ] where = £" 2 00 (/x ~ 
(1/V?) 1/2 ) so d/i = -2r 3 ^ =► (d£/£ 2 ) ~ ~(f/2)dA* ~ -(dp/2y/jZ). Since A x = 0(/i JV //i JV +( 1 / 2 ))^ = 
0{^ {1/2) dn, (l 2 = 0{n 1/2 )diJL (with lead term -(1/2)) we obtain < * >~< T >= 0(1) and 



< A >= O(l). Thus (2.10), (2.11), (C) generate all conservation laws simultaneously with < T 3 > 
(resp. < Xj >) giving rise to f)i (resp. O2). It is then proved that all of the modulational equations 
are determined via the equation 

d T (li = 12<9 x f2 2 ( 2.23) 

where the Riemann surface is thought of as depending on X, T through the points Xj (A, T) . In 
particular if the first 2N + 1 averaged conservation laws are satisfied then so are all higher averaged 
conservation laws. These equations can also be written directly in terms of the Xj as Riemann 
invariants via drXj — SjdxXj for j = 0, 1, • • • , 2N where Sj is a computable characteristic speed (cf. 
d2.32 )-( fj~33 )). Thus we have displayed the prototypical model for the Whitham or modulational 



equations. 

Another way of looking at some of this goes as follows. We will consider surfaces defined via 
i?(A) = ]^[ 2s+1 (A — Ai). For convenience take the branch points Aj real with Ai < • • • < A 2g+ i < 00. 
This corresponds to spectral bands [Ai, A2], • • • , [A 29 +i, 00) and gaps (A2, A3), • • • , (Aa ff , A 2ff+ i) with 
the Ai cycles around the gaps (i.e. ai ~ (A 2 i, A2»+i), i = 1, • • • ,<?). The notation is equivalent to 
what preceeds with a shift of index. For this kind of situation one usually defines the period matrix 
via iBjk — § B u)j and sets 

/ dw.j = 2n5 jk (j, k = 1, • • • , g); duj = ]T (2.24) 
JA k l V H {A) 

(cf. || ^| |l^l)- The Bi cycles can be drawn from a common vertex [Pq say) passing through the 
gaps (A 2 i, A 2 i+i). One chooses now e.g. p — J dp and Q = J dil in the form 

MA) = /*(A) = /^;P = A» + X>A"; (2.25) 
fl(A) = j cffl(A) = J ^—j=3^-dA; = J2 bj A3 ~ j l - -3 ^ A, 



1 
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with normalizations f£ 2%+1 dp(A) — J^' +1 dfl(A) = 0; i = 1, • • • , g. We note that one is thinking 

here of (*) ip = exp[ipx + iClt] ■ theta functions (cf. ( |2.l| ) with e.g. p(A) = —i(logip) x ; fi(A) = 
—i(logip)t- Recall that the notation < , > x simply means x-averaging (or ergodic averaging) and 
(logip) x ~< (logip) x > x ^= here since e.g. {logip) x is not bounded. Observe that (*) applies 
to any finite zone quasi-periodic situation. The KdV equation here arises from Lip = Aip, L = 
—d 2 + q, dt4> = Aip, A — 4<9 3 — 6qd — 3q x and there is no need to put this in a more canonical form 
since this material is only illustrative. 

In this context one has also the Kruskal integrals Iq , • • • , lig which arise via a generating function 

p(A) = -i< {logiP) x > x = ^ + E (2 ^A)2,+i ( 2 - 26 ) 

where I s =< P s > x = P~ s (s = 0, 1, • • •) with -i{logip) x = \ r k + Y^[P s /{2yfX) 2s+1 ]. Similarly 

-z(I^ = - i ^ = 4(VA) 3 + V ^ (2.27) 

and one knows dtP s = d x £l s since (♦) [(logtp) x ]t = [(Jogip)t] x - The expansions are standard (cf. 
ft2[ |l^, [lj], [H| ) . Now consider a "weakly deformed" soliton lattice of the form 

— oo<ni<---<n m <oo j.fc j 

with the Aj (i = 1, • • • , 2g-\- 1) (or equivalently the parameters u 1 — If, i = 0, • • • , 2g) slowly varying 
functions of x, t (e.g. u l = u l (X, T), X = ex, T = et, i = 0, 1, • • • , 2m). Now one wants to obtain a 
version of d]23| ) directly via (4). Thus insert the slow variables in (4fc) and average, using edx or 
£<9t in the external derivatives, to obtain 



d T {logiP) x =dx(logiP)t (2.29) 

or 8tp(A.) = dxE(K). Then from ( |2.25| ) differentiating in A one gets 

d T dp = d x dQ (2.30) 

Now (recall dtP s = d x fl s ) expanding (2.30) in powers of (VA) -1 one obtains the slow modulation 
equations in the form 

d T u s =d x n s (s = 0,---,2g) (2.31) 
where £l s is a function of the u % (0 < i < 2g). This leads to equations 

d T Ai =v i {h l ,---,k 2 g + i)d x k i (i = l,..-,2ff + l) (2.32) 

for the branch points A& as Riemann invariants. The characteristic "velocities" have the form 



dfl 
dp 



A=A 4 P( A *) 



To see this simply multiply ( [2.30 ) by (A — A^) 3 / 2 and pass to limits as A — > Aj 
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2.4 Extension to KP 



Given some knowledge of symmetries as sketched in |Tc|| for example one is tempted to rush now 
to an immediate attempt at generalizing the preceeding results to finite zone KP via the following 
facts. The KP flows can be written as d n u = K n (u) where the K n are symmetries satisfying (in the 
notation of |oJ) the linearized KP equation d 3 (3 = (l/4)<9 3 ^ + Zd{uj3) + (3/A)d~ 1 d^/3 = K'[(3}. The 
conserved densities or gradients 7 satisfy the adjoint linearized KP equation $37 = (l/4)9 3 7+3uc?7+ 
(3/4)9 _1 9|7. Then, replacing the square eigenfunctions by ipip* one nas e -g- i ) '4 > * = So° s n^~ n 
where s n ~ a 7„. Further d n u = K n+ i = ds n+ i = dRes L n — dS7I n+ i where V/ ~ Sf/Su 
(Res L n — nH^ l _ 1 is generally used in the multipotential theory). We are working here in a single 
potential theory where all potentials in L = d + Ui+\d~ l are expressed in terms of 1*2 = u 
via operators with d and One uses here the Poisson bracket {/, g) = J J (5 f / 5u)d(5g / 8u)dxdy 
(Gardner bracket). Let us retrace the argument from [M and see what applies for KP. Thus one 
has s n+ i ~ 7n+i ~ V/ n+ i as conserved gradients satisfying the adjoint linear KP equation (**) 
dtj = (l/4)<9 3 7 + 3udj + (3/4)d~ 1 dyj. The non local term <9 _1 here could conceivably change some 



of the analysis. We need hrst a substitute for ( 2.11 ) or else a direct way of rewriting the adjoint 
equation as perhaps dt[A] + d x [B] + d y [C] — 0. To get such a formula one thinks of differentiating 
the adjoint linearized KP equation to get 

% [dj] + d[id 3 7 + Zudj] + d v [~d y7 ] = (2.34) 

This removes all of the nonlocal terms and one doesn't have to deal with <9[(3/4)<9 -1 c^7] for example; 
however it leads to a redundant situation. Let us try instead 

dth] + d[^d 2 j + Sd^udj] + dyi-dyd- 1 ^] = (2.35) 

One expects that < <9 2 7 >= =< dyd^ 1 ^ > and we take C = 3d~ 1 (uodj) so that, averaging as in 
KdV, one obtains 

d T < 7 > +d x < t >= (2.36) 
which conceivably might be useful. This is discussed later. 



3 AVERAGING WITH 

Let us look now at |5l| but in the spirit of |4(J . We will expand upon this with some modifications 
in order to obtain a visibly rigorous procedure (cf. also fujf ). Thus consider KP in the form ([2~5|): 
3a 2 u yy + d x (Aut — 6uu x + u xxx ) = via compatibility [d y — L, dt — A] = where L = a~ 1 (d 2 — u) 
and A = d 3 — (3/2)ud + w (a 2 = 1 is used in (40| which we follow for convenience but the procedure 
should work in general with minor modifications - note d means d x and u — > —u in the development 
of (|2.l[ )). We have then (d y — L)ip = with (dt — A)ip — and for the adjoint or dual wave function 
ip* one writes in [|l[ tp*L — —d y ip* with ip*A = d t ip* where ip*(fd j ) = (—d) j (ip*f). We modify the 
formulas used in |[l| (and [^0)) in taking 

i> = e px+Ev+nt -(l>(Ux + Vy + Wt,P); ip* = e- px ~ E y- nt ■ <t>*(-Ux -Vy -Wt,P) (3.1) 



where p = p(P), E = E(P), = £l(P), etc. (cf. (2.2) - (2.4)), which isolate the quantities needed 
in averaging. The arguments to follow are essentially the same for this choice of notation, or that 
in j4(J] or |5lj ]. Now one sees immediately that 

(4>*L)^ = ^*L^j + + d 2 x (i>*L 2 ^) + • • ■ (3.2) 
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where e.g. L r = ((— l) r /rl)(d r L / d(d) r ) . In particular 



r xx i> = - 2d(r^ x ) + 2 w>»; (3.3) 
-r xxx ^ = ffex - + 3[5 2 (v*^) - 3{r^ xx )\ 

This means (L* = L, iP*A ~ A*tp* = -d 3 tj)* + (3/2)d(tp*u) + wtfj*) 

(ip*L)ip = [{d 2 - u)ip*]ip = ^* xx i) - uip*ip = tp*Lip - 28(ip*ip x ) + d 2 (ip*i/j) (3.4) 
which implies L 1 = —2d and L 2 = 1. Next 

= -VLx^ + |(V£« + i>*Ux)i> + wip*ip = VA^VO + (3.5) 

- 3<9 2 )?a] + sa 2 ^*^) - a 3 (^» 

so that A 1 = — 3d 2 + (3/2)u, A 2 = 39, and A 3 = —1. Wc think of a general Riemann surface 
E g . Here one picks holomorphic differentials duk as before and quasi-momenta, quasi-energies, etc. 
via dp r~j dfli, dE ~ dft2, dfl ~ • ■ ■ where A ~ fc, p — Jp o dCl x , etc. Normalize the dflk so 



that 5ft J A dfife = = 5ff J B cK7 fe (cf. remarks after (2^)); then U, V, W, ■ ■ ■ are real 2g period 
vectors, and one has BA functions tp(x,y,t,P) as in (3.1). As before we look for approximations 
based on u (xU + yV + tW\I) = u (9j,I k ). For averaging, 6j ~ xUj + yV 3 + tWj, 1 < j < 2g, 
with period 2tt in the 9j seems natural (but note 9j, g +j ~ Uj, etc.). Then again by ergodicity 

< 4> >x— ^TOx,_>oo(l/2L) J^ L 4>dx becomes < 4> >= (V^ 71 ") 23 / '"/ 4>d 29 and one notes that 

< d x 4> >— automatically for <j> bounded. In one thinks of (f>(xU + ■ ■ •) with <fi x — ^ Ui(dcf>/d9i) 
and J • • • J {d(j>/d9 i )d 2a 9 = 0. 



Now for averaging we think of w ~ u {\S\I) as in with (S, I) ~ (5, J)(X, Y, T), 
dxS — U. dyS — V, and drS = W. We think of expanding about uq with 9 X — > 9 X + e<9x- This 
step will cover both x and X dependence for subsequent averaging. Then look at the compatibility 
condition (Jit) : d t L — d y A + [L, A] ~ 0. As before we will want the term of first order in e upon 
writing e.g. L = Lq + eL x + • • • and A = Aq + eA x + • • • where slow variables appear only in the 
Lq, Aq terms. As indicated in [^lj the term [L,A] — ► {L,A} where {L, A} arises upon replacing d x 
be d x + edx in all the differential expressions and taking the terms of first order in e. However the 
formula in pl| is somewhat unclear so we compute some factors explicitly. In fact, according to p0| , 
one can write now, to make the coefficient of e vanish 

d t L x - d y A 1 + [L , A x ] + [L X ,A ] + F = 0; F - d T L - d Y A - (L^xA - A 1 d x L) (3.6) 

Thus F is the first order term involving derivatives in the slow variables (note the L x dxA and 
A 1 dxL terms have opposite signs from |l| but this seems to be correct). To clarify this let 
us write (d - d/dX) L, = (d + ed) 2 ~ (u +«*! + •■•) = d 2 - u Q + e{2dd - u x ) + 0(e 2 ) and 
A, = (d + ed) 3 - (3/2)(tt + em + ■ ■ ■) ■ (d + ed) + w + ewi + ■ ■ ■ = d 3 - (5/2)u d + w a + e{38 2 d - 
(3/2)in<9 - (3/2)m 9 + Wi) + 0(e 2 ). Write then A e = A Q + A x e + 0{e 2 ) and L e = L + L x e + 0(e 2 ) 
with Ax = 3d 2 d - (3/2)uid - (3/2)u d + wi and L x = 2dd - u x . Note that L 1 = -2d ^ L x and 
A 1 = —3d 2 + (3/2)uq A x but we can write 

A x = -A 1 8 - (3/2)u x d + wx = -A X B + A x ; (3.7) 
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L\ = -L 1 d - mi = -L 1 d + Li 
with L\ = —u\ and A\ = — (3/2)md + w\ as in ^]. Then (Jit) becomes 

d t L - d y A + [L Q , Aq] + e{d t Lx - d y Ax + T L - d Y A + [L , i x ] + [L u Aq}} + 0(e 2 ) 

Now the dtLo — d y Ao + [Lq, Aq] term vanishes and we note that 

[L , Ai] + [Li,Ao] = [L , -A l d] + [-L 1 ^ A ] + [Lq, A{\ + [L u Aq] = 

= [L , Ai] + [Lx, A ] + A l dL Q - L 1 6Aq - L A 1 B + A L x d 

Then one notes that OxLq ~ dxL and OxAq ~ dxA (dx ~ <9) so dropping the terms in ([ 
an inoperative d on the right we obtain (|3.6|). Next one writes, using (|3.2|) 



(3.8) 



(3.9) 



with 



dt^L^) - O y (rA^) = r{Lu - A ly + [Li,A] + [L, A x ]}^ = (3.10) 

= ^(dtLx - d v Ax + [L Q , A x ] + [ii, A ])V + &,(■ ■ •) 

and via ergodicity in x,y, or t flows, averaging of derivatives in x,y, or t gives zero, so from (3.6) 
and (3.10) we obtain the Whitham equations in the form < tp*Fip > = (this represents the first 
order term in e - the slow variables are present in Lq, Aq, tJj, and ip*). In order to spell this out in 
pf| one imagines X, Y, T as a parameter £ and considers L(£), A(£), etc. (in their perturbed form) 
with 

(3.11) 



</,(£) = e p({)*+m)y+n(i)t . suffix + ^(£) y + W(0t\I(0) 



where ijj* = exp(—px — Ey — Vlt)<f)*(~Ux — Vy~Wt[I) (no £ variation - i.e. assume p, E, il, U, V, W, I 
fixed). We recall that one expects Xk = Xk(X,Y,T) etc. so the Riemann surface varies with £. Also 
recall that x, y, t and X, Y, T can be considered as independent variables. Now as above, using ([3~2]), 
we can write 

dM*m) = r{A(o - a^o - s^avco) + d l(- ■ •) ( 3 - 12 ) 



(3.13) 



Note also from ( |3.1l| ) for P fixed (9 ~ x?7 + yV + yW) 

W>(£)k=o = (P* + ^» + + 

+(£7x + T>y + Wt) ■ ip*d e if> + i ■ Vdiip 

where / ~ df/d£. In |l0| one assumes that it is also permitted to vary £ and hold e.g. the Ik 

constant whil e allowing say the P to vary. Now differentiate the left side dt(ip*ip(£,)) of ( 3.12 ) in £ 
and use ( 3.13 ) to obtain 

«(V*V>(0)]l£=o = + W ■ We + (3.14) 

+{(px + Ey + {lt)d t (ip*ip) + (tlx + Vy+ Wt) ■ d t (ct>*cj>e) + I ■ W<£/)} 
Fixing x,y,t and averaging ( |3 - 14 ) in the 9 variables yields 

< % WV(0)]le=o >= n < > +^ < ^Vfl > (3-15) 

Next one di fferen tiates the right side of ( |3. 12 ) in £, averages in the 9 variables for x,y,t fixed, and 
equates to (3.15), to obtain 



(3.15) =< >-< 0f[WAV(O)]le=o >= 



(3.16) 
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< 7p*Aifj > -p < ijfA 1 ^ > -U- < ip*A 1 ip e > 

We note here as in (3.13) - ( |3.14| ) (one can disregard A 1 action on the isolated x term, and assumes 
/ fixed, since only more terms d x {- ■ •) arise whose average vanishes) 



0«WAV)|*=o = d x [ip*A 1 (px + Ey + titty + 

+^P*A 1 (Ux + Vy + Wt) ■ ipg] = p^A 1 ^ + U ■ ip*A 1 ->p s + 
+(px + Ey + fLy)d x (i>*A x i>) + (Ux + Vy + Wt) ■ d x (i>*A x ^ e ) 
and ( 3.16 ) follows. Rewriting ( |3.16[ ) with we obtain (note tj}*ipe = 4>*4>e) 

- < Tp*d Y A4> >= -dyCl < ip*ip > -d Y W < cf>*(t>e > - 

-d Y p < ip*A x ij} > -d Y U- < ifA^o > 

Similarly one uses (cf. ( 3.1 2| ) ) 



(3.17) 



(3.18) 



with £ ~ T to get 



Finally note that 



ww)) = rmo - Dm - ^(v*£V(0) + dli- ■ •) 

< ip*d T Lip >= d T E < i)*i> > +d T V- < cj)*<j)g > + 
+d T p < > +d T U- < tp*L x tp e > 

d v {rA i m)-dt{rL i m) = 



(3.19) 
(3.20) 

(3.21) 



= r[L\A{i) -A)- A\L(0 - L)m) + d x {- ■ ■) 
Using this with £ ~ X one gets then as above 

< ip*(L x d x A - A 1 dxD4' >= dxtt < tp*L x ip > -d x E < ^A 1 ^ > 
+ d x W- < ip*L x ipg > -d x V- < i>*A x ip e > 



This formula ( 3.22|) will be of importance later in deciphering the notation of |5E 
dxS — U, d Y S — V, and drS — W so there are compatibility relations 



(3.22) 
We recall that 

(3.23) 



d Y U = d x V; d T U = d x W; d T V = d Y W 
Now add up ( |3.18| ) and ( |3.20| ), and subtract fl3.22| ) to get 

< ^*[d T L - d Y A - Dd x A + A 1 d x L\^ >=< ^j*F^ >= = (3.24) 

= {d T E - d Y n) < tP*iJj > +(d T p - d x Sl) < fl> > +{d x E - d Y p) < ^*A x tp > 
We observe that if one lets the point P on S vary with £, while holding 9k and Ij fixed, then 

df(yV>(£))l£=o = (xdp + ydE + tdn)(f>*(f> (3.25) 



which, together with (3.12) and (3.19), yields e.g. 

dtt < <f>*4> >= -dp < ip*A x <i/j >; dE < 



>= -dp < ip*L x ij > 



(3.26) 
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(cf. - (p7|)). It follows then from ( pM] ) and |T26| ) that 

= (Cl Y - E T )dp + (p T - fl x )dE + (E x - p Y )dn 
Thus one arrives at a version of the Whitham equations in the form 

p T = fix; Py = E X ] E T = fly 



(3.27) 



(3.28) 



where the last equation establishes compatibility of the first two via pty — Pyt — &xy — Ext — 
dxiS^Y — Et) = 0. We feel that this derivation from pQ ], based on |5lj, is important since it again 
exhibits again the role of square eigenfunctions (now in the form ip*^) m dealing with averaging 
processes. In view of the geometrical nature of such square eigenfunctions (cf. JTo| for example) one 
might look for underlying geometrical objects related to the results of averaging. Another (new) 
direction involves the Cauchy kernels expressed via and their dispersionless limits (cf. p8|). 

Let us write now (T ~ px + Ey + fit and H = xU + yV + tW) with 

< >=< ip*(-2d)ip >= -2 < e- T ^*[pe T (t) + Ue T (j)g] > = 

- 2[p < <\>*<\> >+U < <t>*4> e >] 
Similarly from A 1 = —3d 2 + (3/2)u we have 



< >= -3 < tp* [p 2 e r (f> + 2pUe r (/)g + U 2 e r ^ 



+- < ^*u ^ >= 



-3[P 2 < 



> +2pU < 



> +U 2 < 



9]> 



>] + - < <t>*Uo<t> > 



(3.29) 
(3.30) 



We note also a natural comparison with ( 2.22 ) for the formula 

oo 

< S n > 



< ^* >= ■ 



A" 



(3.31) 



We remark here that such a series is natural from asymptotic expansions but when ip, ip* are written 
in terms of theta functions it requires expansion of the theta functions in 1/A (such expansions are 
documented in p7f , p. 49 for example) . It is now natural to ask whether one can express dp, dE, dil 
from ( |3.26| ) in more detail. Thus note first from ( 3.26 ) 



dE < 
3dp[p 2 < ( 



' >= 2dp[p < (j)*4> > +U < (t>*4> 9 >]; dtt < 4>*4> >= 

3 

> +2pU < <f>*(j)g > +U 2 < (jfcpee >] - 7;dp < <f)*U (f> > 



(3.32) 



which says e.g. 



= 3p 2 dp 
E ~ A 2 - 



dE = 2pdp 

< <i>*<i>e 



2U < 



> 



6pU 



> 



< 



3Z7 



> 



< 



dp; d£l = 

g > 3 < (f)*U Q (f) > 



(3.33) 



dp 



< 4>*(p > < 4>*4> > 2 < 4>*4> > 

Given p ~ A+- • • , E ~ A 2 +- ■ • , Q, ~ A 3 + - ■ • this is reasonable. One would also like to compare ( 2.36| ) 
and the equation d T dp+d x dn+d Y dE = (recall £ = (l/A)d 2 -f+3d- 1 (u d-f) and Q = (3/4)^0^7). 
One expects dp = dX + ■ ■ ■ , dE = 2XdX + • • • , and dft = 3X 2 dX + ■ ■ ■ so a formula like (3.31) has 
an order of magnitude compatibility with a possible asymptotic identification dp ^< 7 > g?A. Note 



14 



here the analogy with ( |2.22[ ) etc. where < * >~< T > with ~ - < T > (dp,/2 y /JI) (the 
weight factor arises from KdV notation here, linked to a hyperelliptic Riemann surface); recall also 
< X > is similarly connected to ti.2- Let us assume the "ansatz" (4»4»): For some choice of standard 
homology basis and local coordinate, dp ^< ipip* > dX for large |A|. We have seen that a variation 
on this is valid for KdV situations and via changes in homology and local coordinate it is reasonable 
in general. In fact it generally true via 

LEMMA 3.1. The ansatz XX is valid for KP situations. 

Proof. We refer here to [|l0| for background and notation (cf. also jl3| for dispersionless genus 
zero situations). We write d = L + o~jL~i which implies dtp/ip = A + J^erjA - - 7 . (this is 
motivated by the proof of Lemma 2.1 in 0). Then (logip)x = A + ano - (^ffVO* = p = 

A + J^T < a j > ^ J - But s ™+! = ~ na n ~ dj a n-j an< i one obtains < s n+ i >= — n < a\ > 

- Ei" 1 < 9 X >= -« < <?l >• Thus > dp = d\- J2T J < a ) > A _J_1 cZA = dA + £^ < s J+1 > 
A - -*' - 1 dA =< tpt/f* > since = s n X~ n , and the normalization is built into this construction. 
QED 



Similarly from (3.26) we obtain for large A 

dfl < > dX; dE < ^*LV > dX (3.34) 



(note this requires < -0*L 1 ^ >= 0(2A) which by ( 3.29| ) is correct while < ^)* A x i\j > should be 



0(3A 2 ) which is also correct). It is clear that dp =< ip*i(j > dX cannot hold globally since < ip*ip > 
should have poles at D + D* and the correct global statement follows from |5l], |58|, namely (dp/ < 
>) = dti where dti is the unique meromorphic differential with a double pole at oo and zeros 
at the poles of ip*tp (i.e. at D + D*); this characterizes dti and one can state 

THEOREM 3.2. With the notations as above dp =< ijjtp* > dti, dE < ^L 1 ^ > dti, 

and dti ~ - < ^*A lr ip > dCl. 



REMARK 3.3. From j2~3^ ) one obtains now 8t < 7 >= — dx < C >. Further since pt = tix 



corresponds to dxdp = dxdti and dxdp ~ &t < 7 > dti (note dti will not depend on slow variables), 
one can formally write 

» ^ r ^ » ^ ^ dTd P ^d^ /q 

- dx < C >~ Ot < 7 >~ — — ~ — - — (3.35) 

suggesting that 

< C > dti ~ -dn ~< > dti => < £ >~< > (3.36) 
This is consistent with ( 2.36] ) since 



9 T < 7 > +dx < t >~ 9 T ( * ) - d x ( ^T) = (3.37) 

via drdp — dxdti. 

In summary we can also state 

THEOREM 3.4. The quantity ipip* is seen to determine the Whitham hierarchy a;nd the 
differentials dp, dE, dti, etc. 

REMARK 3.5. Since D + D* — 2oo ~ K^, ipip* is determined by a section of K% (global point 
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of view) but it is relations based on < ip*L Y ip >, < ip*A x ip >, etc. (based on the Krichever averaging 
process) which reveal the "guts" of ipip* needed for averaging and the expression of differentials. 



4 DISPERSIONLESS THEORY 

4.1 General framework for KP 

We give next a brief sketch of some ideas regarding dispersionless KP (dKP) following mainly 
p2| . |l3| , fl4"| , f4l| , fi"7[ |74| to which we refer for philosophy. We will make various notational adjustments 
as we go along. One can think of fast and slow variables with ex = X and et n = T n so that d n — > 
ed/dT n &ndu(x,t n ) — » u(X,T n ) to obtain from the KP equation (1 / 'A)u xxx +3uu x + (3 / '4)<9 _1 d\u = 
the equation d T u = 3ud x u + (3/A)d- 1 (d 2 u/dT$) when e -> (d^ 1 -> (l/e)^ 1 ). In terms of 
hierarchies the theory can be built around the pair (L. M) in the spirit of (h], [l^, |74|. Thus writing 
(t n ) for (x, f„) (i.e. x ~ ti here) consider 

oo oo oo 

L e = ed + Y,u n+ i(e,T)(edr n - M e — ^ nT n L™~ 1 + ^ v n+1 (e, T)L~ n ~ 1 (4.1) 
l 11 

Here L is the Lax operator L = d + J^T u n +id~ n and M is the Orlov-Schulman operator defined 
via ip\ = Mip. Now one assumes u n (e,T) — U n (T) + 0(e), etc. and set (recall Lip = Xip) 



4> 



, oil 



expljr— A") = acp('-S(T ) A) + 0(l)) ; 



i 6 / V £ 



r = exp ( 1f(T) + O ( - ) ) (4.2) 



We recall that d n L = [B n ,L], B n = L%, d n M = [B n , M], [L,M] — 1, Lip — Xip, d x ip = Mip, 
and ip = t(T- (l/nX n ))exp[J2T T n \ n ]/T(T). Putting in the e and using d n for d/dT n now, with 
P = Sx, one obtains 

00 oo 

A = P + ^C7„ +1 P-"; P= A-^PiA- 1 ; (4.3) 

1 l 

oo oo 

X = nT n \ n - x + ^ V n+ i\- n - x ; d n S = B n (P) => d n P = dB n (P) 
i i 

where d ~ <9x + (dP/dX)dp and M — * M. Note that one assumes also (e, T) = Vj+i (T) + 0(e); 
further for B n = J2o ^nmd m one has B n = Y^o bnm,P m (note also B n = L n + ^i^ 7 "')■ We a 
few additional formulas which are easily obtained (cf. jl^l ); thus, writing {A, B} = OpAdA — dAdpB 
one has 

d n X = {B n , A}; d n M = {B n , M}; {A, M} = 1 (4.4) 

Now we can write S = T n \ n + with d m Sj + i — <r™ 1 , V n+ i — —nS n+ i, and d\S — 

M (af Further 

B n = A" + V d n S j+1 \-3; dS n+1 ~ -P n ~ -^±i ~ (4.5) 
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We sketch next a few formulas from p7] (cf. also (4^, 50 1). First it will be important to rescale 
the T n variables and write t' — nt n , T' n = nT n , d n = nd' n = n(d/dT^). Then 

A™ K 
d' n S = -±; d' n X = {Q n , A} (Q n = ^); (4.6) 

d' n P = dQ n = dQ n + d P Q n dP; d' n Q m - d' m Q n = {Q n , Q m } 

Now think of (P, X, , n > 2 , as basic Hamiltonian variables with P = P {X, T n ) . Then -Q n (P,X, T' n ) 
will serve as a Hamiltonian via 

(recall the classical theory for variables (<?,£>) involves q = dH/dp and p = —dH/dq). The function 
S^A, X, T n ) plays the role of part of a generating function S for the Hamilton-Jacobi theory with 
action angle variables (A, — £) where 

A™ 

PdX + Q n dT' n = -£d\ - K n dT' n + dS; K n = -R n = ; (4.8) 

n 

^T = K = d{Rn = 0; J| = & = -d x Rn = -A- 1 

(note that \' n = ~ d' n A = {Q n , A}). To see how all this fits together we write 
HP f)P HX BP 



This is compatible with (4.7) and Hamiltonians — Q n . Furthermore one wants 

Sx = £ S x = P; d'J = Q n - R n (4.10) 



and from (4.8) one has 

PdX + Q n dT' n = -£d\ + R n dT' n + S x dX + S x d\ + d' n SdT^ (4.11) 

which checks. We note that d' n S = Q n = B n /n and Sx = P by constructions and definitions. 
Consider S=S~Z? A n T^/n. Then S x = S x = P and S' n =S' n - R n = Q n - R n as desired with 
£ = S x = S x - J2T Tn^ 1 - It follows that f ~ M - J2T T n A"" 1 = X + Y.T VJ+iA - *- 1 . If W is 
the gauge operator such that L = WdW~ x one sees easily that 

Mip = W (^>2 kx k d k -^j = ^ G + £ ^feA fc_1 ^J ^ (4-12) 

from which follows that G = WxW^ 1 — > £. This shows that G is a very fundamental object and 
this is encountered in various places in the general theory (cf. jHJ |l3|, ?, (79)). 

REMARK 4.1. We refer here also to |l^, Q for a complete characterization of dKP and the 
solution of the dispersionless Hirota equations. 
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4.2 Dispersonless theory for KdV 

We give a special treatment here for KdV since it is needed in [0 |2(| . Thus following || [n], [b| 



we write 

L 2 = L\ = d 2 + q = d 2 - u {q = -u = 2u 2 ); q t - 6qq x - q xxx = 0; (4.13) 

B = 4<9 3 + 6qd + 3q x ; L 2 = [B, L 2 ]; q = -v 2 - v x ~ v t + 6v 2 v x + v xxx = 

(v satisfies the mKdV equation). KdkV is Galilean invariant (x' = x — QXt, t' = t, u' = u + A) 
and consequently one can consider L + d 2 + q — X = (d + v)(d — v, q — A = — v x — v , v = ip x /ip, 
and —ip xx /tp — q — A or ip xx + qtp = Xip (with u' = u + X ~ q' = q — A). The v equation in ( 4.13| ) 



becomes then vt — d(—6Xv + 2v 3 — v xx ) and for A = — k 2 one expands for Sfc > 0, \k\ — > oo to get 
(••) v ~ ik + ( Vn / i^) n ) ■ The v n are conserved densities and with 2 — A = — v x — v 2 one obtains 



ra-l 



P = -2«i; 2v n+ i = - 2J "n-mWm - v' n ; 2v 2 = -v[ (4-14) 
1 

Next for ip" — uip = —k 2 ip write ib± ~ exp(±ikx) as x — > ±00. Recall also the transmission and 
reflection coefficient formulas (cf. |J) T(k)ip- = i?(fc)V>++'0+( — fc, x) and TV>+ = RLip-+i/j-(—k,x). 
Writing e.g. ip + = exp(ikx + (f>(k,x)) with <fi(k,oo) = one has 4>" + 2ik<j>' + (<fi') 2 — u. Then 
^/ + /ip+ = ik + (/)' — v with q — A = — v x — v 2 . Take then 

00 , 

*' = £p^r; t,~ifc + ^ = i* + X;^^^ = 2»« B (4.15) 

Furthermore one knows (cf. §) 

l °a T = - E JSr; c 2»+i = ^- / fe2 "Mi - i^i 2 )* ( 4 - 16 ) 

(assuming for convenience that there are no bound states). Now for c 22 = Rl/T and c 21 = 1/T 
one has as fc — > —00 (Sfc > 0) the behavior tp + exp(—ikx) — > C22exp(— 2ikx) + C21 — > C21. Hence 
exp{<j)) — > C21 asi-> —00 or —00) = —logT which implies 

y< ™ = ? /^ (4,7, 

Hence J (f>2 m dx = and C2 m +i = — J 4>2m+idx / (2i) 2m+l . The C2„+i are related to Hamiltonians 
H2n+i = a n C2n+i as in 0, and thus the conserved densities w„ ~ <^ n give rise to Hamiltonians 
H n (n odd). There are action angle variables P = klog\T\ and Q = r yarg(Rz J /T) with Poisson 
structure {F,G} ~ J (SF/Su)d(SG /Su)dx (we omit the second Poisson structure here). 

Now look at the dispersionless theory based on k where A ~ {ik) 2 = —k 2 . One obtains for 
P = S x , P 2 + q = -k 2 , and we write V = (1/2)P 2 + p= {l/2)(ik) 2 with q ~ 2p ~ 2m 2 - One has 
dk/dT 2n = {{ik) 2n , k} = and from ifc = P(l + qP" 2 ) 1 / 2 we obtain 



1 



ik = P 1 + [ 2 ) q m P~ 2m (4.18) 
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(cf. ( |4.3| ) with ii2 = q/2). The flow equations become then 

dL+iP = SCan+ii = {Q 2n +i,ik} (4.19) 

Note here some rescaling is needed since we want (d 2 + q) 3 / 2 d 3 + (3/2)qd + (S/4)q x = B 3 instead 
of our previous B 3 ~ 4<9 3 + 6qd + 3q x . Thus we want Q 3 = (1/3)P 3 + (\/2)qP to fit the earlier 
notation. The Gelfand-Dickey resolvant coefficients are defined via R s (u) — (l/2)Pes(<9 2 — u)^ 1 / 2 ) 
and in the dispersionless picture R s (u) — ► (l/2)r a _x(— u/2) (cf. fl3fl ) where 



Res 



, ,2^+(i/2) _ ( n+ (1/2) \ n+1 _ (n + 1/2) ■■■(1/2 +1 . 
( ~ fc) n+1 J 9 - (n+l)! ' 

2d g r„ = (2n + l)r„_! (4.20) 



The inversion formula corresponding to (4.3) is P = ik — Y)^° Pj(ik) 3 and one can write 



2 2 

^i-flC + (?) = <9 2 ™+l(- fc2 ); ^2n+l9 = 2n j = ^TT'V™^ = 9X^-1 (4-21) 

Note for example r = g/2, rj = 3g 2 /8, r 2 = 5<? 3 /16, • • • and d' T q = qxro = (l/2)qqx (scaling is 



needed in ( 4.13 ) here for comparison). Some further calculation gives for P = ik — J2T Pn{ik) n 

I />oo 
Pn ~ -v n ~ ~] C 2n+1 = (-1)" +1 / P 2n+1 (X)dX (4.22) 

The development above actually gives a connection between inverse scattering and the dKdV theory 
(cf. 0, [l| [b§ for more on this). 

5 WHITHAM EQUATIONS AND SYMPLECTIC FORMS 

5.1 General relations 

In |5^, |5(| [5?], [5^| t here is a development of Whitham theory and symplectic forms in a general 
context (cf . also |15|, |l6| , [l?], |22|, |3l], [32], [33| ) . The idea is to have a universal Whitham hierarchy 
based on a space M g x of smooth algebraic curves £ s with N punctures P a , a — 1,---,N, and 
local coordinates k~ 1 (P) (here P is a point on S g ). For each such datum and a set of g points 
7i j * " )7s € S ff in general position (or equivalently a point in the Jacobian J(S g )) standard algebro- 
geometric constructions give a quasiperiodic solution of some integrable nonlinear PDE. For KP one 
has N = 1 and for Toda N = 2; we concentrate here on N — 1. Note that P is used to represent a 
point on S ff (corresponding to P = Sx in dispersionless theory), with dp = dfli and p(P) = J P dp. 
From now on we will be scrupulous in distinguishing dQ t and Qi — J P dQi. In fl53[ |5(| one uses 
k^ 1 as a local coordinate at 00 and we will think here of k\ ~ k ~ A where A corresponds to the 
"spectral" variable for KP (which in dKP becomes an action variable). Further p(k, T) is specified as 
fli(k, T) so k and P both represent points on S ff . Thus one must be careful in interperting P = Sx 



from dKP when working on E g . In fact the S of (2/7) satisfies Sx — U as indicated before (3.6) 
and we will see below that there will be an action term S(p) = Y^T'i^iip) with diS — fi*. Hence 
for Oi = p one has dxS = p and this plays the role of P = Sx (from dispersionless theory) on T, g . 
Thus, for simplicity, in what follows P G £ s , k ~ Pi ~ 00 with no other punctures, and p ~ f2i. 
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Now in [|3[ |6| one relates the Whitham equations (44) diQj — djili (di ~ d/dTi), where the 
fii are expressed in suitable variables as indicated below, to zero curvature equations 

8 a ^b - d B 0,A + {Qa,&b} = 0; {f,g} = ,fxg P - f P gx (5.1) 



(note the order of terms in {f,g} is different from ( |2.3| )). Here the A, B correspond to (1, i), 1 ~ P\, 
and i ~ f2j where one thinks of k = k(p, T) or p = p(k, T) with T\ = X, Q,a = &a(p, T), and 



/P~k i 
cMi = k l + Oikr 1 ); k*(p) = W ^P S (5-2) 



The analogy here is to 0, of (4.5) for example where 

oo jj, n 

B n = X n + = X n - V ^i A -" 1 = V b nm P m (5.3) 
1 o 

Note that b n o ^ is possible however and recall that some normalization (£. dOj = for example) 
will hold; in any case the development is parallel. Now if E is an arbitrary function of (p, T) one 
can regard (5.1) as compatibility conditions for 

d A E = {E,n A } (5.4) 

(easy exercise, using the Jacobi identity) and this places one in the context of what are called alge- 
braic solutions (when there is a global solution of ( |5.4| ) . Here global means that E is a meromorphic 
solution of (|J) such that {E(p, T),k(p, T)} = so there exists f(E) with k(p,T) = f(E(p,T)) 
(note then {E, k} becomes Exf'E p — E p f'Ex — 0). Generally one will stipulate that dE be a 
normalized meromorphic differential of the second kind. When d p E ^ one can writ e p = p(E, T) 
and d A f(p, T) = d A F(E, T) + (dF/dE)d A E (F(E, T) = f(p, T)). Then for such an E (O) becomes 



d A n B (E,T)=d B Q A {E,T) (5.5) 
corresponding to (44). Indeed we note that (Q,' A ~ (dClA/dE), E' ~ (dE/dp)) 

(5.1) - d A n B + Q' B E A - d B n A - n' A E B + {dx^A + ^ a Ex)^ b E' - (5.6) 

-n' A E'(d x n B + n' B E x ) = a A n B - d B n A + 
+n' B {E A + E'd x n A ) - n' A {E B + E'd x n B ) = o 

But ([[J) implies for example that E A = E x d p VL A {p, T)-E'd x ^ A {p, T) = E x n' A E'-E'(d x Q A {E, T)- 
n' A E x ) = -E'dx^A E A + E'd x ^ A {E,T) = 0. Similarly E B + E'd x CL B (E,T) = and (fj) 
implies ( |5 . 5[ ) . Thence, via (|5.5|), in coordinates E, T we can introduce a "potential" S(E,T) such 
that 

Cl A (E,T) = d A S(E,T) (5.7) 

and define 

uj = J2n A dT A (5.8) 



(note this S is not a priori related to the S of ( |2.7| )). Then for Q — dS/dE and 5 ~ full exterior 
derivative 

5S(E, T) - QdE = d E SdE + d A SdT A - QdE = tt A dT A = w (5.9) 
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(note also that p = dS/dX = Oi). It follows that 6w — ^25£Ia A <1Ta formally. Note that we are 
treating here dE and cITa as independent objects and in this spirit from OaS = Qa and dS = QdE 
one obtains dQ A = d(d A S) = d A {QdE) = (d A Q)dE so that 

(thus d A dE = in this context - i.e. &AdE ^ d{d A E)). In the same spirit from u> — ^A^dTA — 
SS - S E dE = SS - QdE we get 

Suj = S 2 S - SQ AdE = dE A SQ (5.11) 

This is set equal to SE A SQ in |5^, ||(| but this seems incorrect unless other assumptions are made, 
so a justification follows. Thus first for a more complete discussion of related manipulation we refer 
to [M where (S ~ S generically now) 

M = d x S ~ Q = d E S: B n = d n S ~ Q A = d A S; (5.12) 
dS = MdX + B ndT n ~ SS = QdE + ^AdT A 



We note next that, using (5.10) and d A f — 8aF + FeOaE before (|5.5[), 

{Q(p, T), E( P , T)} = dxQd p E - d p Q8 x E = (5.13) 

= (d x Q{E,T) + Q E E X )E P - Q E E P E X = Q X E P = = 1 

dE dp 

(analogous to {X,M} = 1 in Q, where {A, M} = XpMx - ^xMp). Further note from {Q,E} = 
1, 8aE = {E, £Ia}, and the Jacobi identity one has 

{d A Q, E} + {Q, d A E} = => {d A Q, E} + {{E, n A }, Q) = (5.14) 

= {d A Q,E}-{{Q,il A },E} 

since {{Q, E}, Sl A } = 0. Then from {/, E} = one gets f x E p - f p E x = so for f(p, T) = F(E, T) 
one has (F x + F E E x )E p - F E E p E x = F x E p = (F ~ d A Q ~ {Q, &, A } here). Thus up to functions 
of T n , n > 1 there results 

d A Q = {Q,Q A } (5.15) 



which is not quite satisfactory. A better approach follows from 74 1. Thus if we assume Suj — SEASQ 
one can write out Sui = SE A SQ = A dT A to get 

^2{d p n A dp + d B ^ A dT B ) A dT A = (5.16) 

A 

= (E p dp + J2 dcEdTc) A (Q p dp + d D QdT D ) 



The coefficients of dp A dTA and dX A dTA are then 

_ n f) . B?- 

(5.17) 



d p Q A = E p d A Q - Q p d A E; 



d x Q A = d x Ed A Q - d A Ed x Q 

Solving for 8aQ gives ( 5.15j ). Hence if we stipulate that ( 5.15| ) holds then Suj = SE ASQ is permitted 
in ( jjUjj ). 
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5.2 Times, coordinates, and differentials 

Another step needing further consideration is the construction S(p,T) — J p Q(p,T)dE(p,T) or 
S(p,T) = XX ^4 — E^)Ha(p,T). This works for the genus zero situation with &aS = £Ia (but this 
last formula needs proving - cf. |5(| for a "proof when E ~ A" = A" , A ~ k). We remark also that 
for any formal power series Q(p) = ^iP* one can define times Tj via 

ti = -Mes p=00 k~ l {p)Q{p)dE{p) (5.18) 

Thus for Q ~ (dS/dE) with S{p) = J2T T ^ ~ Ei° T * fc * + 0(k- r ) one has 

1 1 °° 

Ti = -Res k=00 k- l dS = -Resu^^jT^- 1 - 1 = T t (5.19) 

l 

For the general situation without additional punctures we will describe matters as in fl6|| following 
II (cf. also H H IH H). Thus take e ' g - 

E = p n + u n _ 2 p n - 2 + ■ ■ ■ + u + Oip- 1 ) or £ ~ k n + Oik' 1 ) (5.20) 

This includes E ~ A™ or E ~ A™ or dE 1 ~ G?f2„ (when normalized). For Whitham equations we will 
want normalization such as § A dE — but a larger theory can also be envisioned as in ||. Let us 
think of %1§ A . dSl n = = 3t§ B . d£l n here. Similarly take 

dQ = d(k m + 0{k- 1 )) (5.21) 

so dQ ~ dO m is appropriate when normalized. If the dE and dQ are not normalized there are 
g free parameters for each, corresponding to the holomorphic differentials which can be added 
without affecting asymptotic behavior. Look at the configuration space M g (n,m) — {T, g ,Pi ~ 
oo, [k] n , [k] m , E, Q}. Note that an n-jct of coordinates X)o^"' a i ^ s an equivalence class via k' = 
k + 0(fc - ™ -1 ) and we will specify e.g. Q = X)-oo w ith no logarithmic term; the corresponding 
jet [k] m is then determined by m + 1 coefficients Co, • • • , c m plus g parameters for the holomorphic 
differentials which could be added without modifying the asymptotics. Hence in M, g (n,m) we 
will have 3g — 3 + 1 parameters (for the moduli space of Riemann surfaces with genus g with one 
puncture), plus n + m + 2 + 2g parameters for E, Q\ hence there are 5g + n + m parameters involved 
in specifying M g {n, m). This number will be matched by 5g + n + m "times" 

T 3 = -Re Soo {k-iQdE), (j = 1, • • • ,» + m); t m .e = <f dE; (5.22) 

J J Ai 

T~B it E = f dE; t Ai .q = f dQ; r Bi , Q = f dQ; aj = f QdE 

JBi J Ai JBi J Ai 

where i = 1, • ■ • , g. 

Now let T> be the open set in A4 g (n, m) where the zero divisors of dE and dQ do not intersect (i.e. 
the sets {7; dE(j) — 0} and {7; dQ(y) — 0} do not intersect). It is proved in |5j| that near each 
point in V the hg + n + m functions Tk, ta^e, T Bi,E, TAj.Qj T Bi,Q, o-i have linearly independent 
differentials, and thus define a local holomorphic coordinate system. Further the joint level sets of 
these functions (omitting the ai) define a smooth g-dimensional foliation of U, independent of the 
choices made in defining the coordinates themselves. Then Ai g = Ai g (n, m) can be taken as a base 
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space for two fibrations M 9 and TV; Af g has fiber 5 9 (S) ~ J(E) = Jacobian variety (via the Abel 
map (71, • • • , 7 9 ) — * J2i=i Ip d^i) an d N nas fiber S (all over a point (£, P, [k] n , [k] m E, Q) £ M g ). 
We consider only leaves M. of the foliation of T> indicated above and look at the fibration JV or 
AT 9 over the base A4. First, although E and Q are multivalued on the universal fibration their 
differentials are well defined on Af. The idea here is that E and Q are well defined near Pi and their 
analytic continuations by different paths can only change by multiples of their residues or periods 
along closed cycles. But on a leaf of the foliation the ambiguities remain constant and disappear 
upon differentiation. Hence one has differentials SE and SQ on the fibrations which reduce to dE 
and dQ acting on vectors tangent to the fiber. One can trivialize the fibration Af with the variables 
a\, ■ ■ ■ , a g along the leaf M. and e.g. E along the fiber. Then dQ coincides with (dQ / dE)dE where 
dE is viewed as one of the elements of the basis of one forms for Af and the full differential is 
8Q = dQ + J2i(&Q I dctijdcLi = dQ + 6 E Q (also one can take here the constant terms in E, Q to be 
zero). This framework, with QdE = pdE ~ Xsw ~ Seiberg-Witten (SW) differential for example 
leads to a symplcctic form u>m — ]C da-i A duji for Af 9 which can be written in various interesting 
ways (cf. (l6l [fij |H and see also g |, |§ |J, H M> M> H' H M> S 0' the constructions 



appear in Section 5.2. 

In terms of Whitham theory we have differentials dilj = dfa + 0(k^ 1 )) normalized as before 
which are coupled to the Tj plus holomorphic differentials dujj (j — l,---,g) with <f> A dujj = &y 
coupled to the cij. In addition there are differentials dD,E.i and dflQ.i, holomorphic on S ff except for 
the Aj cycles, where they have jumps 

dn+ ti - dn E l = s i:i dE; dn+ z - dn Q t = s l3 dQ (5.23) 

with normalization conditions 



dQ Etj = <b dn Q .j = (5.24) 

Ai ' J Ai 

These are associated to tb u e and tb { .q and we note that for the Whitham theory the times ta^e 
and TAi,Q vanish by normalization of dE and dQ. At this point the most elegant construction of an 
action differential is in [Q for Seiberg-Witten (SW) theory. The construction there is for a Toda 
theory on a hyperelliptic curve with differentials dfl± m but the arguments should be adaptable to 
the finite zone KP framework (see jl6| for a condensed version of this argument in Qj). Indeed the 
moduli in ]44| are treated abstractly (they could be Casimirs for example) and no branch points are 
involved. Then for the complete theory to apply it is only required that the number K of moduli 
hk correspond to the genus g. Evidently this does not happen generally and in the context above 
one would look at level sets of the times Tj, t B it Ei T Bi,Q, T Ai,E, and TA^q, leaving only g variable 
"times" <2j. In the SW theory many situations (e.g. SU(n) susy Yang-Mills) involve constructions 
of curves via flat moduli hk with K — g. In such situations in one constructs by elegant and 
subtle arguments 

g 

dS(h k ,T) = J2 a ^ + £ T " dn « ( 5 - 25 ) 



such that ai = § A dS with 



ddS ddS , . 

dil n ; — — = du>i (5.26) 



dT n 8a. 

There are then dual variables af — § B dS and a prepotential T such that 



•? = £ (5 - 27) 
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For the introduction of "times" <Zj in BA functions etc. we refer to |38|; this will be spelled out later. 

In the Whitham theory sketched above we have indicated the nice features of an S such that 
Q = 8eS but the construction of S with d n S = fl n was left unfinished. We refer to jyl |6|] and 
the discussion of this in Section 6 for a more rigorous picture but will sketch here a few useful 
constructions from |^3|, |56|, Q for some perspective. For the dispersionless theory there is of course 
no problem (cf. Section 4 and || [l2| |l3[ Q Q). In |^3|. a number of theorems are stated 
(mostly without proof) and we will cite a few of these. Thus working on the subspace D of N — 
E, [k] n , Q, [k] m } where dQ and dE have no common zeros one states in |53| that the map 
D — > T = {Ta}, with Ta defined by ( |5.22| ) (without t Au q or t^e), is nondegenerate, so that 
the Ta determine a coordinate system on D. The corresponding dependence of the curve T, g (T) 
and dE(T) defines a solution of the universal Whitham equations. This interesting point of view 
considers the Whitham hierarchy as a way to define a special system of coordinates on the moduli 
space of curves with jets of local coordinates at one puncture (or more). Next, from [ )56| we have 
the following statement. First define IT where dp(Tl s ) = and set tt s = p(Tl s ) (note dp must vanish 
at 2g points since deg(dp) = 2g — 2 - cf. ]72|). Thus p can be used as a coordinate except at the 
2g points II s and the parameters n s , Uf — § B dp then form a full system of local coordinates on 
M* = fc -1 }. The zero curvature form ( p.l| ) is then equivalent to the compatible system of 

evolution equations 

d A k{p,T) = {k(p,T),n A (p,T)} ; d A Uf = d x Uf; (5.28) 



Uf = £ dn A ; d A n s = d A p{Ii s ) = d x n A (Tl s ) 



'Bi 

One can add here also holomorphic differentials du)i as Hamiltonians with corresponding times 
(where dQ ~ dp here). 

If one works with the subspace N of M* with fixed normalized meromorphic differential dE {E 
as in ( 5.20| )) there are two systems of local coordinates. One is given by itj (i = 0, • • • , n — 2), tt s 



and Uf (so 3g + n — 1 parameters whereas one expects 3g — 3 + 1 for Riemann surfaces with one 
puncture and n + l for a normalized n-jet, giving 3g + n — 1 as desired). A second set of parameters 
is given by 

Uf = j dE; E s = E(q s ) (5.29) 

where dE(q s ) = (s = 1, • • • , D — g). Here dE ~ d(k n + 0(fc -1 )) which indicates a pole of order 

n + 1 (k n - 1 dk z- n+1 (dz/z 2 ) = -dz/z n+1 for z=l/k) so deg{dE) = 2g - 2 => the number of 

zeros is 2g — 2 + n + 1 = 2g + n — 1 . For D = 3g + n — 1 we get D — g = 2g + n— 1 so l<s<D — g 
as in |56| . Then the restriction of the Whitham hierarchy to N is given by the compatible equations 

Let now dHi be a normalized differential defined on Ai such that § A dHi = 0. Then there exists 
a unique differential dSn such that dSn is holomorphic on S g except for the Aj where it has jumps 

dS+{P) - dS H {P) = dHi(P) (P £ A,); <j> dS H = (5.30) 

It is then stated in |5(| that for any solution of the Whitham equations on TV there exists constants 
'- a 



r9 and differentials dHi independent of T, such that this solution is given implicitly via 



— (q s ,T) = 0; dS(p,T) =J2(T A - I%)d£l A (p,T) + dS„ (5.31) 
dp ^ 
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(we have added the d in front of S and £Ia)- This is not entirely clear. Then it is stated that these 
relations imply dS = QdE where Q(p) is holomorphic on E s away from the Ai and oo with jumps 

Q+(E)-Q-(E) = ^p-(EeA i ) (5.32) 

on the Ai. This seems formally OK and, in a situation with constant jumps dHi(P) — TQ.idE(P), 
dQ is said to be single valued on S 9 . This seems to say Q + (E) — Q~(E) = Tq^ which corresponds 



to § B dQ — —Tq i by (5.35) below (note one defines the Tq^ differently in the following and we 



refer to Remark 5.1 for some clarification). 

Thus consider now N — dQ, dE} with dQ holomorphic away from Pi ~ oo. Take coordinates 
7r s , Uf, Ui, and the coefficients b\j in Q = Y]^° b\jp^ + 0(p^ 1 ) and use times Tj, aj, Tq^ = — § B dE 

and Ts,i — § B , dQ. The differentials dQ.E.% and dflQj coupled to T E ,i and Tq^ are uniquely defined 
via analyticity plus jumps on Ai of the form 

dfl E l - dfl E i = SijdE; dQ+ . - dQ Q . = SijdQ; (5.33) 



dn E , t = f dn Qti =0 (j = 1, • • • , g) 

As indicated earlier one can consider these times as coordinates on N so that E g = T, g (T), dQ ~ 
dQ(T) 7 and dE = dE(T). Then one asserts in || that for (*•) dS{E,T) = Q(E 1 T)dE one has 
8aS — Ha- In fact one writes 

oo g 

dS = *^2 Tid^i + ( a idui + T E .jdQ,E, 3 + TQ.jdVtq.j) (5.34) 
o 1 

(we have added some d's and the Tqj terms). The "proof is marginal at best and looks at 
8qj only. Consider dQ.jdS(E,T). From (♦•) it follows that dQjdS is holomorphic except on 
Aj. On different sides of Aj the coordinates are E~ and E + = E~ — Tqj so (for E~ E Aj) 
take derivatives in Q(E~ — Tqj) — Q(E~) — T Ei j (cf. ( |5.35 ) applied to E and Ql) to obtain 



d Q,jQ + ~ d Q,jQ~ = {dQ/dE) (i.e. dT EJ /dT Q ^ = = d Qlj Q+ - d Q ^Q~ - dQ/dE) and this implies 
dQjdS(E, T) = dilQj. As a corollary one concludes that E(p, T) and Q(p, T) satisfy the Whitham 
eequations 8aE = {E, fi^} and 8aQ — {Q, Qa} plus the classical string equation {Q, E} = 1. 

REMARK 5.1. We make here a few remarks about differentials with a view toward clarifying 
some matters above. We recall first (following Q) that an abelian diferential is a meromorphic 
one form given locally via u> — f{z)dz where / is meromorphic in its domain. Evidently dw = 
and one defines a primitive via f2(_P) = Jp w on an y simply connected domain. This is called 
an abelian integral and one must be careful in their definition over the whole surface. Thus let 
(Aj,Bj) be a canonical homology basis and define Aj — J A dfl with Bj — J B dtt. Given a cycle 

7 = ^^rriiBi + riiAi) one has J dfl = ^^rriiBi + riiAi). Abelian integrals of the third kind have 
logarithmic singularities and it is necessary to supplement the above cyclic periods by additional 
polar periods in order to specify precisely the multivaluedness. Thus if dfl has poles at points Pj one 
writes Cj — dfl, j = 1, ■ • • ,n, where jj is a cycle homologous to zero containing Pj. Evidently 
Cj = 2iriRes(dQ, Pj) and necessarily J2 S p c j ~ (where sp ~ singular points). One writes S° for 
the surface obtained by removing all A and B cycles. Letting At, Aj and Bj, BJ be the left and 
right edges of the appropriate cuts, S° has boundary OTP = + Bj ~ Aj — BJ). Any abelian 
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integral f2(P) of the first or second kind (i.e. holomorphic or meromorphic with residues equal to 
zero at all singular points) is single valued on S°. At the boundary it is required that 

fi(p)| A+ - n(p)\ A - = -B r , n(P)\ B+ - n(p)\ B - = A, (5.35) 

3 3 3 3 

This allows continuation of Q as a single valued function on the universal cover of E. To determine 
a single valued branch of a third kind abelian integral one must draw additional cuts Tj between P 
and Pj, where Pj are the singular points. Then the polar period formulas are supplemented with 

Q(P)\ T+ - n(P)| r - = Cj (5.36) 

3 3 

The Riemann bilinear relations yield 

(i) J2( A ' k Bk - A k B' k ) = 0; (it) 52(A! k B k - A k B' k ) = (n "i )!dz n-i ^WUo 5 

{iii) ^(ifePfc - A k B k ) = 2TTiJ2 c 3 n '( p j) = Z^i^Cj / dfl' (5.37) 

1 1 1 "^ p o 

where in (i) f2, J7' are integrals of the first kind, in (ii) W is of first kind and fl is of the second kind 
with a single pole at Po with local parameter chosen so that dfi = \(z — Zo)~ n + 0(l))dz, n > 1, 
and in (iii) f2' is of first kind and f2 is of third kind with no more than logarithmic singularities at 
Pj, j = 1, • • • ,n and Cj — Res(dQ, Pj) (the integration contours in (iii) are chosen to be distinct from 
the basic A or B cycles). As a consequence one knows that any abelian differential of the second or 
third kind with zero a-periods (b-periods) or with all purely imaginary (purely real) cyclic periods 
is uniquely defined by its principal parts at the singular points. The basic holomorphic differentials 
dujj are normalized by the condition J A duj k = 6j k and for hyperelliptic situations one can write 

duj =2_,c jk (1 < j < g); c jk = 2ni(A- 1 ) jk ; A jk = / (5.38) 

i A* JA k M 

Here the matrix Bj k = J B du> k is the period matrix and due to the placements of 2iri in the 
normalization one has 5RP < 0. 

REMARK 5.2. In summary, for dKP, S = &nT n is consistent with d n S — B n and 
d\S = M. Evidently d n S = B n and for completeness we show that d\S = M (cf. |ll9| ). Thus from 
the equations before ( ^5| ) 

oo oo oo 

S = J2 T n X n + Y, T "Y, 9 nS j+ i\- j (5.39) 
1 11 

and one recalls that 2T = ^J°T n 9 n F. Hence 

oo oo oo 

d x S = J2 nTnX' 1 - 1 -J2 T »Y. jdnSj+iX^- 1 = (5.40) 
1 11 

OO OO OO 
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But the last series is x ~ j ~ X [ d i S™=i T nPi - Fj] = J2T F j x ~ j ^ which implies that d x S = M. 
For Riemann surfaces with dS as in (5.34) we can write OacIS = cIQa with 



dS — a,i so that this 



is commensurate with ( (5.25 ) for example where Tea = Tqa = 0. Note here tb^e = § B . dE ~ — T ( 



l Q,i 

and TB t ,Q = § B . dQ ~ Tg^ and since these r's are not used much in |5S| ] they should probably be 
renamed to correspond to their natural function here (i.e. interchange the subscripts E and Q). 
In any event it appears that dual variables af to the action like "times" ai could be defined in 
Whitham theory via 



of = J2 Tn f dQn + J2 a i B » +J2 TA f dnA 

1 J B i 1 1 Bj 



(5.41) 



5.3 Symplectic geometry 

We sketch here the construction of Q for the symplectic form u>m following [0. Thus go back 
to A4g(n,m)=_VE g , P% ~ oo, [fc] n , [k] m , E, Q} defined at the beginning of Section 5.2, with times 
defined via ( 5. 22] ). Define T>, Af, A4, and Af 9 as before and work on Af. We recall that although 



E and Q are multivalued on the universal fibration their differentials are well defined on Af. E and 
Q are well defined near P\ and their analytic continuations by different paths can only change by 
multiples of their residues or periods along closed cycles. On a leaf of the foliation the ambiguities 
remain constant and disappear upon differentiation. Hence one has differentials SE and SQ on the 
fibrations which reduce to dE and dQ acting on vectors tangent to the fiber. One can trivialize the 
fibration TV with the variables ai, • • • , a g along the leaf A4 and e.g. E along the fiber. Then dQ 
coincides with (dQ / dE)dE where dE ~ SE is viewed as one of the elements of the basis of one forms 
for Af and the full differential for Q is SQ = dQ + £?(<9<2 y 'da^ddi = dQ + S E Q. 

Now if one considers the full differential S(QdE) on Af it is readily seen that it is well defined 
despite the multivaluedness of Q. In fact the partial derivatives d ai (QdE) along the base AA are 
holomorphic since the singular parts of the differentials as well as the ambiguities are all fixed. 
In particular (d/da,i)(QdE) — duji where du>i is a basis of normalize d hol omorphic differentials 
dujj = Sij with § B duij = bij . To see this note that it is implicit in ( [5.22 ) since by definition of 



'At 



the ai, (dai/daj) = 5^ — § A (d(QdE) / daj) which implies (d(QdE) / ddj) — dujj. Formally then one 
defines u>m on Af 9 via 



UJM 



= <M £ Q(7i)dE(7i) - £ *Q(7i) A dEira) = £ da, A du t (5.42) 



The first expression seems formally reasonable on Af 9 and the last appears to be a calculation of 
the form uim — 8(QdE) — J^ida-i A [d(QdE)/da,i] — J^iddi A du>i. Now go to KP for illustration 
and background. One can work with a (nonspecial) divisor (71, • • • ,j g ) giving rise to quasiperiodic 
functions of t = (t n ) 1 < n < 00, of the form Ui_ n 1 < i < n — 2, 2 < n < 00, which arise 
as solutions of an integrable hierarchy. The BA function is defined as a meromorphic function 
away from P with simple poles at the 7^ (1 < i < g) and an essential singularity at P of the 
form ip(t, z) = exp tn,z~ n ) (l + £i W 2 *) • There is a Lax operator (in general one for each 
puncture) L n = d n + J2o~ 2 u i,nd l with (d/dt n — L n ) i/j(t, z) = (d — d x , x — t\). Thus there is a 
map {S,Pz, 71, • • • ,7 g } — * {ui tfl (t)}. An explicit form for the BA function for KP is given in ]5l| 
(recall 3? § c dft n = for any cycle C here and cf. || for various forms of B A functions - thus the 
form of tp in |58| appears different from (|]^) but it must be equivalent). Similarly the dual BA 
function if)* is defined as before (and denoted by ^ in |pq|). We note that for the Toda lattice one 
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takes N — 2 punctures and a common theme here to some other work is d(Xsw) ~ m a suitable 
sense (cf. Section 5.2). 

An element in N 9 {n, 1) gives rise to a datum in a space N 9 — {wj n (i)|o } via a mapping 

S : (E,P, [z] n ,E, Q,7i,---,7fl) ~> (S, P, z, 71, ' ' ' , 7 3 ) { Wi,n(*)|™T 2 }- Take a real leaf M (i.e. 
dt<f c ,dE = dl<f> c dQ = for all cycles C on T and note this is automatic for dQ ~ dp ~ df2i and 
dE ~ <iQ„ with real normalization). We also take t% ~ a;. One wants now to express lum i n terms of 
forms on the space of functions {u^ n {€)\. First the u^ n {£) can be written in terms of the asymptotic 



z 



BA coefficients £j and one knows that (d x ip/ip) = z +Y^T hsZ s and (d x ip*/ip*) = —z~ 1 + J2T 
(cf. Q and Lemma 3.1, where e.g. = L + ^i° which implies e.g. (dtp/ip) = X + Y^T <rjA 

for A ~ z~ r ). In any event the first n—l coefficients h s = a\ or h* are differential polynomials in the 
Ui t n (initial data £ a (t)\ x=0 — 4>s(t2, • • •) determine £ s for s < n — 1). Writing if s =< h s > one gets 
now p = z^ 1 + ^j* 3 H s z s (recall also from Lemma 3.1 that < s„ +1 >= — n < <j\ >= —nH n ). We 
note that one uses < >^ and < > xy averaging at various places in |58| but generically this should 
correspond to ergodic averaging. 

A result in |M now asserts that for J\f 9 the Jacobian bundle over a real leaf M. of the moduli 
space M. g {n, 1) the symplectic form lo^ can be written as 

UM = Res P <S ^*^^ > dp = 2 < SJ S A j* S*h* n _ s > (5.43) 

Here h s , h* are differential polynomials as above, the differential forms 8h s and 5*h* are defined via 

Sh s = J27=o $ u i,n J2i(-dY {dh s I 'du^l) and, in a standard manner, S*h* s = Yh=o $ u i% XX^s/^ifn)' 
while 5J S is a linear combination of the Sui n , (0 < i < n — 2). Another (essentially more compli- 



cated) form is also given in |58| for the last term in ( 5.43 ) and we will give some evidence for its 
validity without actually proving it (see below). In any event this is a fascinating result but the 
proof in [js8| requires some embellishment which we extract here from ||l6| with some sign changes. 
First one must come to terms with an expression 



\dp J \ < > 



SE = 6p[==-)-[ - - , = (5.44) 



Sp=[^-)dE- (<^* 6L ^>\ ( dp 



dE ) \ tp*ip > J \dE 



where SL n — J^o 2 ^ u i,nd l (see below) and then it will follow from (5.42) (with Sp ~ 5Q and 
SE ~ dE) that 



g 



1 



< 1p*1p > 



"M = >J I I (7s) A dp{ ls ) (5.45) 



The formula (5.44) is asserted to come from |53[ but to see this one has to interpret Sui -n as 



arising from ed x Ui -n (X,T) in the first order term. Indeed we can look at (3.22) with the last 
two terms absent on the leaves of our foliation and written generically for L n ~ £l n as (•♦) < 

r{LidxL k - L\d x L n )^ >= d x Q k < ^*L X J> > -d x n n < ^*L\$ >. Then for L n = L 3 ~ il 3 
and Lk — L\ ~ p one has L\ = d, 8L\ ~ d x L\ — 0, L\ = —1 and — < ip*d x L 3 ip >= d x p < 
^*L\i\) > +9x^3 < ip*?P >■ Written in terms of Sp ~ ed x p, SL 3 ~ ed x L 3 , etc. one obtains 
(5^3 = -Sp(< i>*L\i\) > / < ijj*t/j >) - (< ip*SL 3 ip > j < >) while from Theorem 3.2 (with 
L\ — A 1 ) we have (< ip*L\ip > / < ip*tp >) = -(dCl 3 /dp), or generically (dfl n /dp) = -(< ^L^tp > 
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/ < ip*ip >), and (5.44) follows. The formula (•♦) can also be used for a more general theorem 
in [ p8[ . Thus ( |5.4i: ) holds and writing 7 S (t) for the zeros of ip away from P (corresponding by 
Riemann-Roch to the fixed poles 7 S for t = 0) one can pick t\, ■ ■ ■ ,t g (generically) as times for which 
the flows 7 S (t) are independent and use them as coordinates on 5 s (E). These can be transferred 
to the system of coordinates /(71), ■ • ■ , f{l g )) for / an Abelian integral on E via (d/dti)f(j(t)) = 
i?es 7 ( t )[(9 / 'dti)ip(t, z) /ip(t, z)]df . In the present case, writing (Sti/i/ip) = WVO^j and using 

Res ls — —Resp one obtains ujm — Resp [(< ip*SL n tp > / < ip*ip >) A {Stip/^P)\ t=0 \ dp. It is not 
unnatural to see here an apparent change in the number of parameters. 



Now note that ip*5L n %l) = ip* J^o 2 5u t ^ n d l ip 
h s z s from which (d^ip/ip) = z~ J [l + c 
>= [-z~ 2 + 0(l)]dz and V*V> = 1 + Ei° s « z 



En — 


.ZP+ 1 ] 



' &Ui_ ri %l)*d l i\) and we recall dip/ip 



+ 



J ( c jp = c jp(h s , d l h s )). Then since dp/ < 
for small z, we have 



■4i*SL n i{j ' 



dp 



J2g 2 Su in tfj*d l ipdp 



(5.46) 



-5> 



z p+1 1 («- a +o(l))(l+> 'a n « n )«fa 



1 



00 



1 



1 



-n— 1+s 



where 5J S is by construction a linear combination of the 6ui n , (0 < i < n — 2). Averaging now 
.46) and using (5.44) one obtains < SJ S >— for s = 1, • • • , n — 1. In addition it is claimed 
that 8H S = for s = 1, • • • , n — 1, since the mean values H s =< h s > are fixed along M. 
(this is not entirely clear). Hence lom = Ylx < 5Jj+n > A dtj. Further from jIo|, |l3| one has 




(d 3 r/r 



djij}*(8L n il>) _ djjS ^ 



This leads to 



dp 



(5.47) 



— n — 1+s 



(fa 



which implies 



n-l 



(5.48) 



Hence ujm = Resp(< <5 t V>* A > / < "0*^ >)^P h olds a nd an argument is suggested in 

|p8[ to extend <5 f to a full 5 so that the first formula in ( 5.43 ) is prov ed (we have made some 
adjustments in minus signs from jl6[ Q). For the second formula in ( [5.4S ) one looks at the relation 

Slogip* = 6 (c(U, i>2)+ f* Q d x logip*J= ( c s(U, i>2)+ f* o h* s z s dx} . This implies 



Res P { <(Slo 9 r)^SL^>^ dp = UM= ^ < SJs A |* 



5*h^_ s dx > 



(5.49) 



which is the second identity in (|5.43| ). In |58| one observes that by definitions SJ S does not contain 
variations of derivatives of ui and then claims that 5J S = —n(dh s /du)Su ~ —n6h s . This seems 
unclear but we will write out some relevant formulas below. Finally we note that there is a little 
interplay between S ~ edx and X = ex in the last two lines which we do not elaborate upon here (cf. 
|l9f for further details of calculations, etc.). The examples in [BS| have some curious features but 
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in any event one should keep in mind the conditions SH S = when selecting functions. One notes 
that of course the symplectic structure involving 6J S , S*h*, etc. (i.e. the 8ui n ) does not include the 
time dynamics so KdV ~ L 2 — d 2 + u for example and Bousinesq ~ L 3 = d 3 + ud + v. Adaptions 
and applications of the above fibration framework to SU(n) 1 N = 2 susy YM theory are also given 
inf|. 

In order to compare the 8J S and Sh a we make a few calculations here. Thus 5J S is given via 
as a linear combination of the Sui„ with complicated coefficients depending on h s , d l h a and the 
coefficients s n in ip*ip (note the s n can be expressed via & and £* while the term dp/ < ip*tf) > can be 
written (z~ 2 +J2^ PmZ m )dz where the p m are independent of the u, n ). We recall here (cf. & [To|, |65| ) 
for -0 = exp(£)[l + J2T ij zj \ and ^* = ex P{~0[ l + J2T witn L = d + J2T u i+i®~ 1 i u 2 ~ «) 

one has u 2 = u 3 = — £ 2 + e tc- Further from dip — z^ 1 ^ + exp(^) J2£i zl and (9ip/ip) — 
z^ 1 + h s z s we get 

00 00 00 n— 1 

(1 + £ Ci^xE = E ^ fe & = h m + E ( 5 - 5 °) 

111 1 

so e.g. hi = £i, £ 2 = Zihi + h 2 , C3 = h 3 + h 2 £i + &hi, etc. Hence hi = -u 2 , h 2 — C 2 - £ihi = 
£i£i ~~ u 3 + £i u 2 = — U3, etc. In the present situation one is interested in L™ ~ L n so we note that 
L\ = d 2 + 2u 2 , L\ = d 3 + 3u 2 d + (6u' 2 + 3u 3 ), etc. (we recall also L n = d n + £[p 2 u^d 1 ). Thus 
for n — 3 we have U03 = 3«3 + 6u 2 , U13 = 3u 2 , etc. Calculations based on this seem to agree with 
the formulas in J58|, namely, for fixed n 

hi = -; h 2 = u n _ 2 ; • • • (5.51) 

n n n 

(formulas for h* are also given in Thus for n = 3 one has hi = — (l/3)iii3 = — u 2 and 

h 2 = (2/3)it' 13 — (l/3)wo3 = — u 3 + 2u' 2 — u 2 = —u 3 , etc. As for the Q s we write e.g. 

^ = [z -i + STh s z°f + J2h> s z s = (5-52) 
™ 1 1 



z- 2 ' 



2 [1 + 2^ h s z s+1 + h sh k z s+k+2 + h 2 s z 2s+2 + Y, Kz s+2 \ 

s^k 

Then writing this as z~ 2 [l + c 2p z p+1 ] one has (•♦) c 2 \ = 2hi, c 22 — 2h 2 + h\, c 23 = 2h 3 + 
h\ + h' 2 , ■ ■ -. Now look at ( 5.46[ ) for n — 3 so we have to consider 5uo3ip and Sui^dip, yielding for 
S = tp*SL 3 ip the formula 

3 = [5u 03 + Sursiz- 1 + £ h s z s )](l + ^ s m z m ) (5.53) 
Putting in the (dp/ < ip*ip >) term we obtain for n = 3 

oo 

(5.46) = Y 8J s z s - A dz = [5u 13 z- 3 + Su 03 z- 2 } x (5.54) 
1 

oo oo oo 

X [1 + J2 S rnZ m ] [1 + ]T h s Z S+1 ] [1+J2 PkZ k+2 ]dz = 
110 

= [5ui 3 z~~ + 5u 03 z~ 2 ]{l + siz + (s 2 + hi + po)z 2 + • • -]dz 
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Consequently SJ± = Suis — 3£it2 = — 3Sh\, 5J 2 = siSuis + Suqs, etc. Note s% = £1 + £* = so 
8J2 = Suqs = (5(3^3 + 6u' 2 ) = —35h 2 + 66u' 2 (note also h 2 ~ — 3^03 mod u 2 ). We see that in some 
"generic" manner the relation S J n ~ —n5h n seems to be holding (modulo terms involving derivatives 
of the Ui n ). Hence provisionally one could say that there is some evidence for a formula of the form 



n-l 

lum r~j —n 



n-l „ x 

^2<Sh s A / 6*h* n _ s > (5.55) 
1 ^ 



as in [p8f , where Sh s only deals with the appropriate potentials and not their derivatives. However 



the expression (5.43) in terms of SJ S is simpler and more direct so we treat this as the basic such 
formula. 



6 REMARKS ON BA FUNCTIONS 



We go here to |jq] for a construction of BA functions with certain a, variables inserted explicitly 
which connect to the a, "times" of SW theory via ieaj = aj. This leads to expressions for S and 
F (action and prepotential) in two different forms from which some comparisons can be made, of 
use in various directions (cf. jl7|, 0, p8|). The construction in J68| is for the Toda theory (with 
two punctures on the Riemann surface) but, in keeping with the spirit of this paper, we will display 
matters for the KP situation only. The idea is to insert new times oti in the BA function (2.1) (recall 
Bij — § B duji, U = (Uj) = [l/2m) § B dfi 1 , etc.) via for example 



ip = exp 



rP a 

I xdft 1 + ydVL 2 + tdfl 3 + i V] onduji 

JPn 



(6.1) 



6 (A(P) +xU + yV- 



Q(A(P) + zo) 

This is an alternate version of the BA function in |68[ |, modulo t n — (cf. H) for relations) or 
one can adapt the formulas of Q to write (recall A(P) — (Jp o dujj) and we set for convenience 
Ojfe = J B d£lj where dflj ~ d£l J ) 



ip = exp 



dW ' +0 J '(P ( 




dujj + uij(Po) 



(6.2) 



e {A(P) + Erfe/ 2 ™)(^) + 1 E? *j(B jk ) + z ) Q(A(P a 



zo) 



e {A(Poo) + ETi-tj/^K^k) + * E? + Zo ) &(A(p) + z ) 

(note one can write here also J P dttj ~ f p dQj + Qj(Po)) and A(P) = (Jp dujj) + A(Poo), where 
flj ~ A J — Ei 2 \<lm,j l'm)z m ~ 1 )- We have attached a factor of (l/2m) to the in order to utilize 

(S.18) below. This kind of formula (|S.2|) provides a normalization as z — > (A — > 00) since, for any 



Pq, the theta function quotient tends to one. Next one constructs in an ad hoc manner the algebraic 
form of ip near 00 via vertex operators as (A ~ z^ 1 ) 



ip = exp 



T{t-[\- l \,g) 
r(t,a) 



(6.3) 
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(here [A x ] ~ (l/n\ n )). Note that the on j duJi terms do not contribute to the essential singularity 
at oo. 



Now the general action and prepotential theory associated with ( |6.1| ) (cf. [T^, [l9| |44|, 68 ) 

leads to 

,„ , v^rr, ,^ ddS , <9<iS I( _ 

do = y cijdLUj + y L n d\l n \ — — = du>i\ „„ = ctsi,, 



1 1 
If we consider functions F(a,T) related to dS via 



<9a, 



2?ri 



9T„ 



-ReSocZ n dS 



(6.4) 



then, given the standard class of solutions of the Whitham hierarchy satisfying 



OF 



<9F 



there results 



2F = J2^~ f dS ~ £ T n Re Soo z- n dS 
2m J B . 



(6.5) 



(6.6) 



(6.7) 



(connections to pa for example will be spelled out below). Writing now, in the notation of [B8j 



^i°°>n 2: " 1 l dz with cifi n = [— nz 



dujj 



expanded to give 



Eoo 
1 In 



]dz, and using (3.4), (2.8) can be 



^ 9 g oo oo 

2F = — ^2 B jka ] a k + 2^ a ^ a jkF k + £ qkiT k Ti 



(6.8) 



j,k=i 



k,l=l 



Thus the expression ( |6.8| ) comes from the Riemann surface theory, without explicit reference to 
the BA function, and we consider now (S.2) and (6.3) to which ideas of dKP can be applied to 
introduce the slow variables T k . This means that we will be able to introduce slow variables in two 
different ways and the resulting comparisons will show an equivalence of procedures. In practice 
this will enable one to treat e on the same footing in the Whitham theory and in the dispcrsionless 



theory, and subsequent analysis will verify a construction in |17| . Thus from (6.2) and ( |6.3| ) one 
obtains an expression for r of the form (t% — x, t%=y, t% = t, • • •) 



r(t,a) = exp[F(a,t)]<d ^A(Poo) + ^fo/27rz)(% fe ) + aj{B jk ) + z j 
where k = 1 , • • • , g and 

Y x - y x 9 x " x 

F ( a , f ) = ^ £ qkitkk ~ — £ B jk aja k + i^Ja.,- £ Vjktk + £ dkh 



fcj=i 



(6.9) 



(6.10) 



j,k=i 



(see also |l5| for a similar form - recall here A{P) = (Jp dujj) and P Q =/= P m is required - note 
that we have inserted a factor of l/2m with the (Qj k )). Putting in the slow variables T k — et k and 
a k = iea k one will find that the quadratic part of F in T is exactly F(a,T)/e 2 . One recalls here 
that r will have the form (4.2) so that (F(a,T) /e 2 ) + (J2dkT k /e) is a natural object (modulo a 
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multiplicative factor of lower order arising from the thcta function in ( |6.9| ) , which could be absorbed 
in an exp{A/e) term). To say this another way, (|4.2| ) is the natural form of r based on ( |6.3[ ) and it 
is associated with tp ~ exp[(l/e)S + 0(1)] as in (4.2). Given the Riemann surface background based 
on (|6.2|), one is expressing the 0(l/e) term in r = ea;p[(l/e 2 )P + 0(l/e)] as exp[(l/e) ^fcPfc] times 
a suitable asymptotic expansion of theta functions (for such expansions see ^7j, pp. 40-49). In |6S| ] 
one writes then from (6.9) and (|6.l[) respectively 



/o - 9r (7'£) =e_2 £ e " F(n)(T ' a); 



(6.11) 



T a 



where dS^ - dS in (|j) and ~ P in (g 



Let us spell out some of the details now. Thus one notes in fl6.2| ) that as P — * P>o the theta 
function terms go to one and the asymptotic behavior is mainly determined by the exponent. Now 
consider identifying fl6.2| ) and (6.3) as P — > P^ (P — ► Poc ~ z — > ~ A — > oo) via (Q asy ~ theta 
quotient in (|6.2D) 



(6.12) 



(write Wj(Po) = J^f dwj and note that d(z m ) — mz m x dz = —mX m 1 d\). If we accept an ansatz 
- ( |6~To|) , then e.g. (p^) can be written out as indicated below. We write i J2i a j {Bjk ) + zq = 
z(a) with z (a) + J2T / 27 ™)(%fc) = z(a,t) and AP = P(a,i- [z]) - F(a,t). Now look at AP 
from ( |6TT0| ) as 



LXJ 

AP = i y 

2 ^ 

k,l=l 



<lkl 



z 

tk ~T 



1 x 



k.l=l 



tl ) — tkti 



z l z k 

tk J + tl T 



From ( 6.12 ) we want then to balance 



kl 



S °° z k 00 z k 

- i Y,J2 a w k T~^2 dk T = (6 - 13) 



3=1 fc+1 
9 oo 



3=1 fe=l 



Y. dk '- 



logQ a 



^ gin " ~fft 



(6.14) 



-z 0ff e ^(Poo) - 



1 



.7 



2?ri 



+ z -log e^PoojH- 5) 



Note also that A(P) — A(Poo) ~ (Jp dwfc) ~ — a km(z m /m)- The balance demanded then 
involves at first order a balance of the explicit tk terms with the remaining terms lumped together in 
dk- In particular this seems to require that qu = (which is true since the sum of residues 

of a meromorphic differential is zero - cf. fl6.20 )) so that J2lkitk{ z l /I) — J2tkJ2 1ik(z l /I). Then 
the explicit % terms in AP cancel the term ^2 tj Qmj{z m /m) in ( |6.14| ) and we can determine dk 
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provided the theta function terms have a suitable expansion. 

In fact one can look at the theta function balance by directly comparing 



A = 



D 



e(^(Poo) + 5) 
e(A(p) + s)e(A(p co ) + z ) 



(6.15) 



e(^(Poc) + z)Q{A(P) + zo) 

First we note some formulas from [||, [l5| [l^, [7l| . Using the Riemann bilinear relations one can show 
that if dn n = Or™ + 0{z))dz then 



d\l n - il nj - - — - 1 yO z UJ j \ z=0 - - _ jT<Tj,n- 



(6.16) 



where dwj = fjdz = — Oj m z m 1 dz — — Gj, p+ \z p dz with ujj ~ — Oj m {z m /m). Note 
also that d n_1 //j| z=0 = —<Jjn{n — 1)! so 



6 -_ 27r V -_ 27r V-if 

ifi n+l,7 — u in — i u J 7 

,J n n! J 



(6.17) 



2=0 



as in [|18| for example. Now in j68[ one writes 

df2„ = Cl n j — —2i:iResz~ n du)j = 2muj n (6.18) 

and this is correct. To see this go to the derivation in [jn] where one shows by the Riemann bilinear 
relations that 

<£ d£l n = Cl n j = 2, ( f duij f dfl n — <£ dujj <£ dfl n ] = (6.19) 

Jb 3 \JAi JBi JBi J Ai J 



2niRes(ujjdfl n — 2niRes ^— <jj m ^ • (^—nz ™ 1 — 



QmnZ 



dz = 2man 



Note also at this point that since the sum of residues of a meromorphic differential is zero, one can 
write 

Res Q m d£l n — Qnm Qthti 



Res 



(6.20) 



showing that q mn = q nm . Further we remark that Res d(u}jdQ n ) = = Res dujjQ n + Res u>jdQ n 
Now to deal with ( 6.1 5| ) one can write 

z 3 (fijfc) 



j 2ni 



+ z 



(6.21) 



Since Sljf. = 2ma^j and — J p duj^ — —uj k = X) (T km{z m /m) we have a cancellation ^ &km(z m /r 
X Ofcj (z-'/j) = for 1 < A; < g. It follows that 



O \AiPoo) - ji^jk) + Sj = G(A(P) + z) 



(6.22) 
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and 

A ejAjP) + go) 

B 9(^(F oo )+z ) 1 j 

This will then lead to dk — dk(a). Indeed ( |3.14| ) becomes now (after cancellation of the tj terms) 

1 °° z k+l a 00 z k °° z k f A\ 

k,j=i 11 1 v / 

Writing «4(P) = -A(Pqo) + (Jp duj) and referring to Q one has for I = ^(Poo) + z 

/ rP \ a 00 m 

logQ ( / dujj) + z - Zog6(z) = - V VVZo^e^ ;(^)cr im — + • • • (6.25) 
V / ,- =1 1 m 

Putting this in ( |6.24| ) yields dfc = dk(a) (cf. [§, |l8| for another version of this). 

For perspective however let us make now a few background observations. First we refer first to Jl^ | 
where it is proved that F mn — F nm in B n = A" — J2T(Fnm/ m )^~ m (the F mn being treated as alge- 
braic symbols with two indices). Since near the point at infinity we have f2„ ~ A™ — ^^°(q mn /m)A^ m 
the same sort of proof by residues is suggested (F mn = ~Res\[B n dX m }) but we recall that B n = A" 
so there is an underlying A for all B n which makes the proof possible. Here this is not a priori 
evident. For example (•♦•) p = A — (^i/ j)^ ^ corresponds to P = A + Pj + iX~ j in Q 
with Pj+i = Fij/j (i.e. Hj ~ —Fij) and the "inverse" is A = P + YT U n +iP~ n (arising from 
a Lax operator L via dKP) The corresponding inverse for (•♦•) then characterizes A in terms of 
p but one does not expect Q„ ~ A" . The matter is somewhat subtle. Indeed the BA function is 
defined from the Riemann surface via d£l n , ckvj, and normalizations. It then produces a unique 
asymptotic expansion at 00 which characterizes tp near 00 in terms of A and hence must characterize 
the dQ n and duij asymptotically. Moreover the normalizations must be built into these expan- 
sions since they were used in determining ip. But the formal algebraic determination of B n via A" 
is a consequence of relating the dVt n to operators L n = L™ as in |p^| . Applying this reasoning 
it is suggested that (•♦•) is valid as a consequence of the BA function linking the differentials 
and the asymptotic expansions. Note that it is not stated that qij = Fji when Fij ~ didjF, 
and this point will be clarified below. When slow variables T n are introduced the moduli hk will 
be functions of the T n (h^ ~ branch points for hyperelliptic curves). This means that a priori 
the normalized differentials dtt n and duij can depend on the T n , along with Bij, etc. The con- 
struction of dS = y~) T n dCl„ + ^2a,jdujj such that d n dS — dVL n and d n a,j — is somewhat subtle 
and we will look at this below following [jl6| [4~4|]. Let us also say a few words about the con- 
struction (3.5), where we recall the formulas ( 5.22j ) for example for the times Tj. In the present 



situation we simply think of dS = QdE in which case ( f]ggh becomes Tj = (l/^Res^X^dS. 



Note here from dS — ^T„c?il„ + ^djduij with dfl n = [riA™ -1 + 0(X~ 2 )]dX one obtains immedi- 
ately ReSoo(l/m)\~ m dS — T m . As for the formula d n F = ReSooX n dS = —Res^z^ n dS we recall 
first the form of S in dKP, namely from Section 4.1, S = y_]T„A" + Sj+xX^ where in fact 

S j+1 = -(l/j)djF or simply S j+1 (l/j)Fj for Fj an algebraic symbol (cf. @ and gg|). Then 

in the dKP case dS = (J2 nT n X n ~ l + OjFX'^dX so Res o0 X n dS = d n F ~ F n is immedi- 
ate. The present situation with dS = ^ T n dCl n + Y a-jduij, dCt n = [nX 71 ^ 1 + ^ q mn X~ m ~ 1 ]dX, and 
= <Tj m X~ m ~ x dX gives ResX k dS = Yl^nQkn + Yl a -i°"-ifc and this is better treated indirectly 
as indicated below; note this suggests Ffe m = qk m (cf. [|18|). 

Let us make some comments around the development in B4] (cf. also |lq] ) . For situations arising 
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in Seiberg-Witten (SW) theory it is desired to find dS — ^T n dVl n + ^ajdojj with d n dS = dVL n 
and d n aj = so that (T„, Oj) can be taken as independent variables with moduli hk = hk{T, a) (1 < 
k < K - note generally K > g). Here the moduli automatically depend on the slow variables T n and 
hence a priori so do the differentials. First one notes that differentials dfl n , having the same asymp- 
totic properties as the dQ n (but unnormalized) , can be selected via dCt n — d[X n — Y]]° cfcm(A~V*)] 
with ai n "arbitrary" and T independent. Then we write du>j = — Pj m z m ~ 1 dz — ^2j m \~ m ~ 1 d\ 
where the <7j m depend on h — (hk) (and thence on T). Now a general desideratum in SW theory is 
that ddS I dhk = ^2 Pjkdujj and to achieve this one picks first the dtl n to have this property. Thus the 
ai n become functions of hk (hence of T) and we can specify ddCl n /dhk = J2i °fei t ^ w » (f° r essentially 
arbitrary a r k \ - modulo the production of analytic functions, which seems to allow enough flexibility 
here). Note dd(l n /dhk = [J2i(9cti n dhk)X~ 1 ~ 1 ]dX which can be written as a linear combination of 
the duij for each n via 



da; 



dhi 



= 5> 



(6.26) 



This defines dai„/dhk and hence the ai„ via the crj^ ; a priori any (reasonable) collection of crjy 
should do. To relate the dCl n to d£l n we write then 



dCl n = dfl n + c^duoi] c™ = j> dCl n 



(6.27) 



Now impose the requirement d n dS = d£l n while assuming dS — ^2u m (T)dCl m (which guarantees 
that ddS/dhk — Pjkdojj). This leads to 



d n ds = 



du ™ ,a \ - ddh m 

-^-dil m + u m dh 0„ii:, 



(6.28) 



E 



dfl r , — dCl r . 



} c n idui 



This implies 



— T 



(6.29) 



E 



d n hk ( 



E 



rp m 
J-m&ki 



^d n h 



k&ki 



where a k i = J2r 



T n m 
± ™- a ki- 



The equation (444) J^^nhk^ki 



-c" here represents the Whitham 



dynamics for the moduli hk- We note that the aj m (h) are inflexible but the terms cL are flexible 



they determine the c™ from (3.27) and the o^i as in ( 6.29| ). The Whitham dynamics are then 
"compatible" equations in some sense - given a choice of a r k \ they describe a compatible variation 
of moduli; however they may not be determining for hk unless K = g with ofo an invertible square 
matrix. Now write 



dS 



Trndfl 

I. 



T 



l dujj 



(6.30) 



CLidlOj 



E T " 



dS 



and look at the requirement d n a.i — which would yield as an independent variable with h k — 
hk(T n ,Zi) along with ddS/don = duii. This situation is handled in via an assumption K = 
g leading to determination of the d n hk via (444) d n hk = — ^2 c f&ik ■ Additional calculation 
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yields then (dhk/daj)\T=c = &jk{h)- One can also examine matters however without making any 
assumption on K. Thus e.g. given (6.3C), d n a,j — requires 



(6.31) 



d n dj = d„ (p dS = (p dQ n = 



which is automatically true for this standard normalization of dtt n . Therefore without explicitly 
solving for hk{T n , aj) we can say that this form prevails when the "Whitham equations" (4k4k4t) f° r 
hk have a solution (perhaps not unique). In these circumstances we will have also ddS/dai = dcui 
and 

af = j> dS = j> aiduji + ^ T n dfi n J = ^ ajBy + ^ T n Q, n j (6.32) 

3 

Let us develop this a little further following [|58| . We will simply repeat some of the argument 
there since it is instructive and revealing. We assume the dfl n and dujj are specified as before 
with aTt n = (-nz _n_1 — J2T q mn z m ~ 1 )dz and dujj = — Vj m z m ~ 1 dz. The standard Whitham 
equations based on times T n and aj are (see Section 5) 



d n dQ m — d m dfl n ; d n dujj — 



ddQ n ddui 



j 



3 daj ' dai 
This form automatically leads to a differential dS such that 

d n dS = df2„; — — = du>j 
dan 



dduji 
daj 



(6.33) 



(6.34) 



for which ( 6.33 ) descri bes compatibility conditions. Then one is led (automatically) to co nstruct a 
function F satisfying (5.5). The consistency of this stipulation (|6.5| ) is ensured by ( |6.33 ) and the 
Riemann bilinear relations. For example to prove d n {dF / da^) = d(d n F) / dai using the formulas of 
( |6.5[ ) one can write (fi„ = J dfl n ) 



(6.35) 



-t 

2iri ^ 



d£l n <t> diOi 

A; JBi 



B, 

dui i> dfl n 

Aj JBi 



= —ReSoaUJndLUi = —Res z n diOi 
(cf. (|5~37l ) and flsls|)). Finally one shows that d n d m F — d m d n F via 



-d n Res z m dS — —Res z m dVt n = q mn — 



Res z" 



Qpn 



dz 



(6.36) 



whereas d m Res z n dS — q nm . From the relation ( |6.5| ) we see also that the local behavior of dS can 
be described in terms of F via 



dS 



Y,nT n z- n - 1 -J2d n F z n ~ l 



(6.37) 
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